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ABSTRACT 

Many current and future astronomical surveys will rely on samples of strong gravita- 
tional lens systems to draw conclusions about galaxy mass distributions. We use a new 
strong lensing pipeline (presented in Paper I of this series) to explore selection biases 
that may cause the population of strong lensing systems to differ from the general 
galaxy population. Our focus is on point-source lensing by early-type galaxies with 
two mass components (stellar and dark matter) that have a variety of density profiles 
and shapes motivated by observational and theoretical studies of galaxy properties. 
We seek not only to quantify but also to understand the physics behind selection bi- 
ases related to: galaxy mass, orientation and shape; dark matter profile parameters 
such as inner slope and concentration; and adiabatic contraction. We study how all 
of these properties affect the lensing Einstein radius, total cross-section, quad/double 
ratio, and image separation distribution, with a flexible treatment of magnification 
bias to mimic different survey strategies. We present our results for two families of 
density profiles: cusped and deprojected Sersic models. While wc use fixed lens and 
source redshifts for most of the analysis, we show that the results are applicable to 
other redshift combinations, and we also explore the physics of how our results change 
for very different redshifts. We find significant (factors of several) selection biases with 
mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile 
inner slope and concentration; concentration of the stellar and dark matter depro- 
jected Sersic models. Interestingly, the intrinsic shape of a galaxy does not strongly 
influence its lensing cross-section when we average over viewing angles. Our results 
are an important first step towards understanding how strong lens systems relate to 
the general galaxy population. 

Key words: gravitational lensing - galaxies: photometry - galaxies: structure - 
galaxies: elliptical and lenticular - galaxies: fundamental parameters 



MOTIVATION 



The field of strong lensing is undergoing a period of growth 
that is expected to accelerate in the coming decade. Cur- 
rent sample s with tens of lenses, such a s from CAS- 
TLE ^ (e.g.. lFalco et~alll200lf). CLAS S (e.g.. iBrowne et all 
20031). SPSS (e.g. , Ilnada et all I2008D . and SLACS (e.g. ; 



Bolton et all l2008h . will give way to samples with hun 
dreds or even thousands of strong lensing systems dis- 
covered by Pan-STARRS, LSST, SNAP, and SKA (e.g. 
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CASTLES is a collection of uniform HST observations of mostly 
point-source lenses from several samples with differing selection 
criteria, rather than a single, uniformly-selected survey. 



Fassnacht et al.ll2004|; iKoopmans et all 12004 : iKuhlen et ail 
20041 : iMarshall et al.ll2005l ). As lens samples grow, they will 
be even more useful for constraining the mass distributions 
of the strong lensing galaxy population (for a review of 
strong lensing astrophysics, see lKochanek|[200r3) . 

There are, however, certain complications in using 
strong lensing studies to constrain the physical properties 
of galaxies. The issue we address here is selection bias. Sup- 
pose we treat galaxies as having an intrinsic probability dis- 
tribution in some parameter x (such as the inner slope of 
the density profile), and we want to use observed strong 
lens systems to infer that distribution p(x). If the lensing 
probability itself depends on x, then in general the distribu- 
tion psl(x) for the strong lens galaxies will not reflect the 
true, underlying p(x) for all galaxies. The more diverse the 
galaxy population, and the more the strong lensing cross- 
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section varies across the population, the more important the 
selection biases can be. 

Obviously, one must account for selection biases in or- 
der to draw reliable conclusions from current and future 
strong lens samples. Our purpose here is to use a complete, 
end-to-end, and above all, realistic simulation pipeline for 
strong lensing (presented in Paper I of this series, van de 
Ven, Mandelbaum, & Keeton 2008) to quantify and under- 
stand strong lensing selection biases and assess their impact 
in typical situations. We focus on strong lensing of quasars 
by early-type, central (non-satellite) galaxies at several mass 
scales, using two-component mass profiles (dark matter plus 
stars) that are consistent with existing photometry and 
stacked weak lensing data fr om SDSS (|Kauffmann et all 
120031 ; iMandelbaum et aT]|2006h as well as with iV-body and 
hydrodynamic simulations. 

The issue of selection biases in strong lensing is not a 
new one, and we are building on many previous studies of 
this issue. The main galaxy properties that we study in an 
attempt to understand selection bias are 



• Mass 



ll99ll ; lMoller et all 12007); 



Turner et all Il984l ; iFukugita fc Turner! 



Shape (e.g. , Keeton et alj 1991 
iRozo et al.ll2007n ; 



• Orientation ( "inclination bias" , e 
iKeeton fc Kochanekll 19981; IRozo et al.lT2007h: 



Mailer et al.1l 19971 ; 



Huterer et all 2005; 



• Dark matter inner slope (e.g. . IKeeton fc Madaul l200ll ; 
IWvithe et aill200ll ; iLi fc Ostrikerll2002l ) ; " 

• Dark matter concentration (e.g., iKuhlen et al.l [20041 ; 
iFedeli et alj|2007l ). 

In most cases, previous studies addressed the issue of se- 
lection bias using models that were simplified in some way. 
Simplifications often included testing effects of dark mat- 
ter slope and concentration using pure (generalised) NFW 
profiles without a baryonic component; or testing effects of 
density profile or shape using single-compone nt mass mod- 
els. A n otably different app roach was taken bv l Hilbert et all 
i|2007l ). lM511er et all (|2007l ), and lHilbert et alj (|2008l ). who 
used the Millennium simulation with a semi-analytic model 
of galaxy formation to simulate strong lensing. That ap- 
proach naturally yields a realistic distribution of halo and 
galaxy properties (to the extent that the true distribution 
of galaxy properties and their relation to halo properties 
is encoded in the semi-analytic model), at least for a fixed 
cosmological model, and it will undoubtedly be useful for 
modeling the strong lensing population in large, future sur- 
veys. We elect, however, to take a more controlled approach 
here, using realistic but discrete values of the density profile 
parameters, so that we can disentangle the different types of 
lensing selection biases, and understand the physics of each 
one, before recombining them to determine the net effects. 

Our basic approach is to generate realistic galaxy mod- 
els with various values of the relevant parameters, and to 
investigate how the lensing cross-section o(x) depends on 
model parameters x (where x might include mass, shape, 
concentration, and other parameters, some of which may be 
correlated). The reason for focusing on the cross-section is 
that it represents the weighting factor that transforms the 
intrinsic joint parameter distribution p(x) into the distribu- 
tion psl(x) among observed strong lens systems. Note that 
we focus on selection biases related to physical effects. There 



may also be observational selection biases (e.g., iKochanekl 
Il99lf ). but they are specific to a given survey and are not 
something that we can address in a general way. This fact is 
one reason for our focus on point-source lensing: extended 
source lensing is even more sensitive to observational selec- 
tion effects than quasar lensing. For example, the effect of a 
finite aperture is no longer straightforward for an extended 
source, and detection of an extended strong lens system de- 
pends on factors such as the size of the point-spread function 
which varies spatially and temporally in any given survey. 
Fortunately, with large numbers of point-source lens sys- 
tems anticipated in fu ture surveys [e.g., ~ 10 5 with SKA 
l|Koopmans et al.ll2004f )] this investigation of physical selec- 
tion biases in point-source lens systems should be highly 
useful. We emphasize that our purpose is only to determine 
the mapping function a(x) due to physical differences be- 
tween galaxies, and not to predict the observed distribution 
of properties psl(x) in some particular lensing survey, for 
which knowledge of the intrinsic parameter distribution p(x) 
and of any observational selection effects are both necessary. 

In this paper, we present a systematic investigation of 
physical strong lensing selection biases that both unifies and 
extends previous work. We begin in Section [2] with a brief 
review of the simulation pipeline developed in Paper I, in- 
cluding the density profiles and shapes of our galaxy models, 
our computational lensing methods, and the basic lensing 
properties of our galaxies. In Section [3] we investigate se- 
lection biases related to galaxy orientation and shape. In 
Section 21 we consider selection biases related to the inner 
slope of the dark matter component of the density profile 
(for cusped models) . In Section [5] we allow the dark matter 
concentration to vary as well. In Section [U] we turn to de- 
projected Sersic density profiles and examine selection biases 
with both the stellar and dark matter parameters. Section]?] 
includes an exploration of the ranges of lens and source red- 
shifts for which our conclusions are applicable. Finally, in 
Section [51 we summarise our main conclusions and discuss 
their implications for past and future strong lensing analy- 



2 SIMULATIONS 

Here we summarise the strong lensing simulation pipeline 
that was developed in Paper I. For more details regarding 
the choice of galaxy models, and the lensing calculations 
(including validation and convergence tests), we refer the 
reader to that paper. 



2.1 Notation and conventions 

The pipeline assumes fiducial lens and source redshifts of 
zl = 0.3 and zs = 2.0, respectively, although we explore a 
range of lens and source redshifts in Section [7] We compute 
distances using a flat ACDM cosmology with Q m = 0.27 and 
h = 0.72. With these choices, the angular scale at the lens 
redshift is 1" = 4.36 kpc and the lensing critical density is 
E c = 2389 M pc -2 . 

We use x, y, and z to denote intrinsic, three-dimensional 
(3d) coordinates, and x' and y for the projected, two- 
dimensional (2d) coordinates. Similarly, r is the intrinsic 



radius (r 



+ y' 2 + z ) and R' is the projected radius 
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(R' 2 = x' 2 + y' 2 ). For non-spherical galaxies, we use a, b, and 
c for the major, intermediate, and minor semi-axis lengths 
of the intrinsic, 3d density profiles; and we use a' and b' for 
the major and minor semi-axis lengths of the projected, 2d 
surface densities. Subscripts "dm" and "★" are used to indi- 
cate whether a quantity describes the dark matter (DM) or 
stellar component of the galaxy model. Any gas component 
that has not cooled to form stars — which is expected to 
be sub-dominant r elative to the stellar component for early 
type galaxies (e.g., iRead fc Trenth am 2005; Mor ganti et al.l 
I2006T ) — is implicitly included in the component that we la- 
bel dark matter. 

All quoted dark matter masses Miso are defined using 
comoving quantities such that the average density within the 
virial radius is 180p. This definition is in principle arbitrary, 
but this choice was used to analyse the weak lensing results 
that we use to estimate masses for our galaxy models. 

2.2 Galaxy models 

Our guiding principle in designing the galaxy models and 
choosing their parameters is that we want them to be real- 
istic examples of galaxies, but not necessarily to reproduce 
the full intrinsic parameter distribution of real galaxies. The 
reason for the latter point is that strong lensing is inherently 
non-linear, so it may be the case that a region of parameter 
space with very few galaxies is actually quite important for 
strong lensing because it has a large cross-section. A classic 
example is brightest cluster galaxies, which are a negligi- 
ble fraction of the galaxy population, yet constitute a non- 
negligible fraction of strong lenses. Consequently, we choose 
parameter values taking into account prior knowledge of the 
intrinsic parameter distributions — where known — and us- 
ing physical intuition for what effects may strongly increase 
the lensing cross-section, such as a steep inner slope of the 
dark matter density. 

2.2.1 Basic setup 

All of the ellipsoidal galaxy models used for this pipeline 
involve two components, one for the stars and one for 
the dark matter. The parameters of stellar components 
are based on photometr y and spectroscopy from the SDSS 
(|Kauffmann et al"1l2003l ). The parameters of the dark mat- 
ter components are derived from JV-body simulations and 
stacke d weak lensing observation s of early-type galaxies in 
SDSS (|Mandelbaum et all 120061 ). We attempt to approxi- 
mately bracket the range of luminosities of observed strong 
lensing systems by using models corresponding to ~ 2L* and 
~ 7 L« early-type, central (non -sat ellite) lens galaxies. (S ee, 
e.g jKeeton fc ZabludoflhoO^I and lMomcheva et al.ll2006l for 
a discussion of environmental effects when the lens galaxy 
is a satellite in a group or cluster of galaxies.) The deriva- 
tion of the resulting mass and length scales is explained in 
Paper I. 

Here we briefly mention that the stellar components 
have masses M to t,* of 1.16 and 5.6 x 10 1 M©, and projected 
half-mass radii -R^,* of 4.07 and 11.7 kpc for the above lower 
("galaxy") and higher ("group") mass scale, respectively. 
The dark matter components have halo masses Miso.dm of 
0.47 and 9.3 x 10 13 Mq, and concentrations Cdm of 8.4 and 



5.6, which implies scale radii r Si dm of 49.7 and 201 kpc, for 
the lower and higher mass scale, respectively. We note that 
a simulation with only two mass models is not sufficient for 
a fully quantitative study of mass bias in strong lensing. 
We have limited the number of mass models because mass 
bias is among the oldest and best-studied selection effects in 
strong lensing, so we focus instead on other effects for which 
the answers are less clear. 

2.2.2 Density profiles 

We use two basic types of density profiles for the simulation 
pipeline. The first are cusped density profiles, which have 
the general form 

p{r) = m ~<{i + m y-^ (1) 

where m is a dimensionless ellipsoidal radius; in the spheri- 
cal case, m — r /r s for some sca le radius r 3 . For the stel- 
lar cusped density, we use the iHernquistl (|l990l ) profile, 
which has (7,71) = (1,4). Fo r the cusped densit y of the 
DM halo, we use the NFW l|Navarro et all 1 19971 ) profile, 
which has (7, n) = (1,3), but generalise it by allowing 7 to 
vary on a grid of values {0.5, 1, 1.5}. Further variation of 
the cusped models includes several simulations with adia- 
batic contraction of the DM halo due to the presence of the 
stellar component , acco rding to the prescr i ption s given in 
iBlumenthal et al.l (Il986l ) and lGnedin et al.l (|2004l ). 

T he seco nd type of profile that we use is a deprojected 

ISersid il 19681) profile, which is well-described using the for- 
mula from Tprugniel fc Simienl (|l997l ), 

p{m) = exp \-b n m 1/n ] , (2) 

where the central (negative) logarithmic slope is related to 
the Sersic index n by p n = 1 — 0.6097/n + 0.05563/n 2 , and 
b n = 2n - 1/3 + 4/(405n ) + 46/(25515n 2 ) to high preci- 
sion ilCiotti fc Bertinl 19991 ). The scale radius is the projected 
half-mass radius, which in case of a constant stellar mass-to- 
light ratio corresponds to the effective radius R e . The Sersic 
indices n * of the stellar comp onent from observational con- 
straints l|Ferrarese et al. 2006) and rtd m of the dark matter 
component from JV-body simulations ( Merritt et al.l 120051 ) 
are 4.29 and 2.96 (respectively) for the lower mass scale, 
and 5.87 and 2.65 for the higher mass scale. We also inves- 
tigate varying both nd m and away from these fiducial 
values within the ranges derived in Paper I. 

Paper I includes a number of figures characterising the 
intrinsic (3d) and projected (2d) density profiles of our 
model galaxies. We highlight the following projected quan- 
tities (since the projected mass distribution is what gov- 
erns lensing): the surface mass density E(-R') (figure 4); the 
negative logarithmic slope "/'{R') of the surface mass den- 
sity (figure 5); the projected enclosed dark matter fraction 
/ dm (i?') (figure 6); and the lensing deflection angle a(R') 
(figure 7). Those quantities are plotted for both the cusped 
and deprojected Sersic density profiles at both mass scales. 

One characteristic of these galaxy models that was 
noted in Paper I, and that will be important in our in- 
terpretation of results in this paper, is that from ~ 1-50 
kpc (cusped models) or from ~ 0.3-30 kpc (deprojected 
Sersic models) , the composite galaxy model is approximately 
isothermal. This result can be most easily seen in Paper I 
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figures 5 and 7. Moreover, it is still nearly true not only 
for our fiducial models but for our models with varying 7d m 
(cusped) or ridm and (Sersic). For the higher mass mod- 
els, the deviations from isothermality are more significant. 

While we will find later on that the deviations from 
isothermality on small scales, well within R e i n , can have 
noticable effects for the lensing properties (e.g., the rela- 
tionship between the Einstein radius and the total unbiased 
cross section differs from the relation a — 7ri? 2 in that ap- 
plies to singular isothermal spheres), we can nonetheless use 
this approximate result to reason about selection biases with 
certain galaxy properties. For example, steeper dark matter 
inner slopes 7d m and higher concentrations Cd m , ridm, or n* 
mean that a larger fraction of the total galaxy mass (stellar 
plus dark) is located within some fixed aperture. When com- 
bining this fact with the approximate isothermality, we can 
envision the profile on those scales as simply being increased 
by a multiplicative factor when we increase 7dm, Cd m , rid m , 
or n t . Consequently, higher values of these parameters in- 
crease the total lensing cross-section, and in the case of 7dm, 
Cd m , and ndm, they also increase the projected dark matter 
fraction f dm in the inner regions. 



Monte Carlo techniques. Specifically, mockcsec places many 
random sources behind the lens, solves the lens equation to 
find all the images, and then tabulates the cross-sections for 
lensing with different numbers of images (2, 3, or 4)0 We 
can add those cross-sections to find the total lensing cross- 
section, otot- We use bootstrap resampling to estimate the 
statistical uncertainties in the cross-sections. 

We account for magnification bias in various ways in- 
tended to mimic different survey strategies. We use the 
GRAVLENS routine mockLF to define an observationally- 
motivated double-power-law luminosity function for the 
source quasar population; and we consider surveys with 
faint, intermediate and high flux limits (ranging from lo- 
cally shallow to locally steep parts of the luminosity func- 
tion: 0.04, 0.4, and 41/*). We compute magnification bias 
using either the total magnification of all images, which is 
appropriate for surveys that do not initially resolve strong 
lens systems, and separately, using the magnification of the 
second-brightest image (for surveys with sufficient resolution 
to resolve strong lens systems from the outset). All told, we 
examine six types of biased cross-sections in addition to the 
unbiased lensing cross-section. 



2.2.3 Shape models 



2.4 Basic results 



We consider seven different models for the shape of the 
galaxy density: spherical, moderately oblate, very oblate, 
moderately prolate, very prolate, triaxial, and a mixed 
model. In all but the mixed model, the axis ratios used for 
the dark matter and stellar components are related in a sim- 
ple way (with the stellar component being slightly rounder 
than the dark matter component). In all cases the axes of the 
two components are intrinsically aligned. The mixed model 
includes an oblate stellar component that has short axis 
aligned with the short axis of a triaxial dark matter halo. 
The chosen axis ratios for all models, and the method of 
sampling in viewing angle to compute accurate orientation- 
averaged lensing cross-sections, are given in Paper I. 

2.3 Strong lensing calculations 

Here we briefly summarise the strong lensing calculations, 
which are described in far greater detail in Paper I. In all 
cases, we consider only point-source (quasar) lensing, and do 
all calculations us ing an update d version of the GRAVLENS 
software package (lKeetonll200Ibl ). 

For the cusped models, we compute the scaled surface 
mass density k = S/E c on a map, and use Fourier trans- 
forms to solve the Poisson equation and find the lens po- 
tential. We then compute the first and second real-space 
derivatives of the potential in Fourier space, and use inverse 
FFTs to transform back to real space. Paper I presents tests 
to determine what map resolution and size are necessary to 
avoid numerical errors in the Fourier analysis. With depro- 
jected Sersic models it is hard for the map-based approach to 
resolve the steep central profile, so we instead compute the 
lensing deflection and magnification analytically (for sp her- 
ical models; also see ICardonel 120041 : iBaltz et all 120091 ). or 
with numerical integrals (for non-spherical models). 

We use the routine calcRein to calculate the Einstein 
radius i? c i n for general density distributions, and the rou- 
tine mockcsec to compute strong lensing cross-sections using 



Here, we present the basic lensing properties of our spher- 
ical cusped galaxy models. Table Q] gives the Einstein radii 
and unbiased lensing cross-sections (a to t = o"2 for spherical 
models) for models at both mass scales, for different values 
of the dark matter inner slope 7d m - 

There are three points to note about these results. First, 
the strong lensing properties are, as expected, a strong 
function of mass. For the fiducial dark matter inner slope 
7dm = 1, the Einstein radius R c i n and lensing cross-section 
increase by factors of 2.6 and 3.1 when we go from the 
lower to the higher mass model; for 7d m = 1-5, the fac- 
tors of 4.3 and 10.7 increase are even l arger. Our results ar e 
broadly consistent with co n clusions by Turner et al.l l|l984h . 
iFukugita fc Turned <|l99ll) . iMoller et all (|2007l ) and from 
many other studies that the population of lens galaxies dif- 
fers from the overall galaxy population in its overall higher 
mean stellar and DM halo mass. The higher mean mass for 
the strong lensing systems is determined by the competition 
between the rising <T to t with mass and the steeply declining 
mass function at high masses. A second point is that the in- 
crease in otot with mass is notably larger for the models with 
7dm = 1-5 than for the models with 7d m = 1- This means 
the strong lensing mass bias depends sensitively on the inner 
slope of the intrinsic density profile. A third point is that, 
even when the mass is fixed, -Rem and er to t are very sensitive 
to 7dm- While this is not surprising, the latter two points 
do serve as a reminder that calling strong lens galaxies a 
"mass-selected" sample of galaxies is an over-simplification. 

2 In each 2-image and 4-image system, there is actually one ad- 
ditional image that lies near the centre of the lens galaxy (so 
in some of the lensing literature such systems are called 3-image 
and 5- image systems instead). The central images are highly de- 
magnified and rarely observed, however, which is why most ob- 
served lenses are considered to be doubles and quads. We do find 
the central images but do not include them in our astrophysical 
analysis. See Paper I for more discussion. 



© 0000 RAS, MNRAS 000, 000-000 



Strong lensing selection biases 5 



Table 1. Summary of the basic lensing results for the spherical shape model. 



Description 




Lower 


mass 






Higher 


mass 






7dm = 0.5 


7dm = 1 


7dm = 1-5 


SIS 


7dm = 0.5 


7dm = 1 


7dm = 1-5 


SIS 


-Rein (kpc) 


2.5 


3.0 


4.5 


2.4 


5.1 


7.6 


19.4 


17.3 


R e ; n (arcsec) 


0.58 


0.68 


1.04 


0.55 


1.18 


1.75 


4.44 


3.97 


""tot (arcsec 2 ), unbiased 


0.42 


0.53 


1.40 


0.95 


0.84 


1.63 


15.0 


49.5 



As a reference point, we also consider a single- 
component singular isothermal sphere (SIS) model which 
has the same mass within the virial radius as our two- 
component galaxy models (for each mass scale). Even if the 
density profile is truncated at the virial radius, for radii well 
inside the virial radius the intrinsic density profile can be 
given as psis(r) = Mi8o/(47rri8or 2 ), and the surface density 
as £(-/?') ~ Mi8o/(4ri8o-R')- To find the Einstein radius, we 
require for the mean surface brightness that E(i? i n ) = E c , 
which gives 



Lower mass scale 



R. 



(SIS) 



2ri 8 ()E c 



oc M 



2/3 



(3) 



where Miso and riso are related by our virial mass definition 
3Mi 80 



180p 



(4) 



The resulting Einstein radii and lensing cross-sections 
(ertot = 7r7iei n ) are also given in Table [1] 

Comparing the familiar SIS model with our more com- 
plex, two-component models yields some notable results. For 
the lower mass scale, the SIS Einstein radius is similar to 
that of our shallow 7d m = 0.5 cusped model, but the SIS 
cross-section lies between those of the steeper 7d m = 1 and 
7dm = 1.5 cusped models. For the higher mass scale, the SIS 

(SIS) 

Einstein radius R cin is only slightly smaller than that for 
the 7dm = 1.5 cusped model, and the SIS cross-section is 
considerably higher than that for any of our cusped models. 
This result can be understood from paper I figure 5, which 
shows that while 7' for our lower mass models is near the 
SIS value of 7' = 1 near i? e i n , it is significantly shallower 
for a fair fraction of the range inside R e iu- The shallower 
inner slope moves the (source-plane) radial caustic inwards, 
decreasing the cross-section to well below the SIS value of 

7Ti?e in . 

One implication is that a simple SIS model cannot 
mimic any particular one of our cusped models. Thus, even 
though galaxies seem to have quasi-isothermal profiles in 
the vicinity of the Einstein radius (see Paper I), a sim- 
ple SIS is still not very realistic as a global lens model. 
The difference between the SIS and more complex two- 
component lens models is particularly notable on mass scales 
correspondin g to groups of galaxies a nd above (also see 
Keetonl Il998l; iPorciani Madaul l200(il ; iKochanek fc Whitel 
200ll ; lKuhlen et al.ll2004F T 



3 ORIENTATION AND SHAPE 

We begin our main science results by examining selection bi- 
ases as a function of galaxy orientation and shape, for fixed 
dark matter inner slope 7dm- In Section [3.11 we investigate 
"orientation bias" for each of the six non-spherical shape 



Very oblate 
(a'/a) dm =l 



Very prolate 
( b '/a) dm =0.6 




0.4 0.6 
(b'/a) to 



0.6 0.8 
(a'/a) 4 



Figure 1. The total lensing cross-section <T tot (top) and 
quad/double ratio (bottom) as a function of projected semi-axis 
lengths for the lower mass cusped model for very oblate (left) and 
very prolate (right) shapes. The lines indicate magnification bias 
modes as follows: black solid = unbiased, red dotted and blue 
dashed = weighted by total and second-brightest magnification 
respectively for 0.04L* limiting magnitude, green long-dashed and 
yellow dot-dashed = the same for 0.4L* limiting magnitude, and 
magenta dot-long dashed and cyan short-long dashed = the same 
for 4L* limiting magnitude. 



models. In Section 13.21 we consider orientation-averaged 
cross-sections for all model shapes, to see if there are sig- 
nificant selection biases associated with the intrinsic shapes 
of galaxies. Finally, in Section [373] we consider whether there 
is any simple, universal conversion between unbiased and bi- 
ased lensing cross-sections. 



3.1 Orientation 

In this subsection, we fix the dark matter inner slope 
to 7dm = 1 and examine how the orientation of non- 
spherical models affects their lensing cross-sections ("ori- 
entation bias"). Our basic physical intuition suggests the 
following trends with orientation for a given shape: 

(i) When viewing the galaxy along the long axis (with 
scale length a), the projected scale radii are minimised, so 
that the most mass is squeezed into the smallest projected 
area. The resulting higher surface density at fixed R' im- 
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Triaxial shape 



Mixed shape 




D.75 O.BO 0.85 0.90 0.95 0.75 0.80 0.85 0.90 0.95 

Figure 2. Contour plots of the total lensing cross-sections CTtot 
(top) and quad/double ratios (bottom) as a function of projected 
semi-axis lengths for the lower mass cusped model for triaxial 
(left) and mixed (right) shapes, where the axis labels indicate the 
shape of the dark matter component. All quantities are unbiased 
cross-sections only. 



plies an increased Einstein radius and higher lensing cross- 
section. 

(ii) The four- image cross-section 04 and the quad/double 
ratio 04/02 increase with the projected ellipticity; they are 
maximised when viewing the halo along the intermediate 
axis, such that the long and short axes are seen in projection 
(i.e., projected axis ratio b' /a! = c/a). The four-image cross- 
section vanishes when the projected surface density has cir- 
cular symmetry, i.e., for an oblate galaxy viewed face-on or 
a prolate galaxy viewed end-on. 

(iii) The three-image cross-sections can only be signifi- 
cant in cases of extreme projected ellipticity, when the tan- 
gential caustic protrudes outside the radial caustic. 

We first consider the very oblate and very prolate 
models. Fig. [T] shows the total cross-section <T to t an d the 
quad / double ratio 04 / 01 for the lower mass scale (results for 
the higher mass scale are similar). The total cross-section is 
indeed higher for oblate models viewed edge-on rather than 
face-on, and for prolate models viewed end-on rather than 
side-on (i.e., viewed with the long axis along the line-of-sight 
in each case). The total cross-section varies with orientation 
by factors of ~ 2 for these strongly non-spherical models; 
in the higher mass model (not shown), the variation is as 
much as a factor of 4 because the flatter DM component is 
more prominent. This finding implies that if strong lensing 
systems seen in reality are drawn from significantly non- 
spherical galaxies, we should observe preferential alignment 
of the long axis with the line-of-sight, rather than random 
orientations. 

The cross-sections for three- and four-image systems be- 
have in accordance with our expectations. As can be seen 
from 04/(72 in the bottom panels of Fig. [T] four-image sys- 
tems are plentiful for orientations that yield a large pro- 
jected ellipticity — for oblate galaxies viewed edge-on, or 
prolate galaxies viewed side-on — but they naturally disap- 
pear as the projected mass density approaches circularity. 
While we have not plotted 03 since in most cases it is zero, 
we see substantial numbers of three-image systems only for 
edge-on oblate galaxies, where the projected matter density 
is so flat that the tangential caustic protrudes beyond the 



radial caustic. Above some critical projected axis ratio, 0-3 
becomes consistent with zero. 

Magnification bias increases the effective cross-sections, 
because lensing magnification brings in systems that would 
have been below the flux limit had they not been lensed. The 
effect is strongest when the flux limit is in the steep part of 
the source luminosity function (the 4L„ limiting luminosity). 
In all cases, the biased cross-section is larger when calculated 
using the total magnification rather than the magnification 
of the second-brightest image, but the trends with viewing 
angle are typically preserved. The 0-4/(72 ratio is increased by 
magnification bias, because four-image lenses tend to have 
higher magnifications than two-image lenses. This ratio in- 
creases with the flux limit, but it is actually higher when 
using the second-brightest magnification rather than the to- 
tal magnification, because many four-image lenses have at 
least two bright images. We discuss further details of the 
biased cross-sections in Section [3.31 

Next, we consider the more complex triaxial and mixed 
shape models shown in Fig. [5] For both mass scales, the total 
cross-sections have maxima when a /a = b/a and b' /a — c/a 
(i.e., the projected semi-axis lengths have their minimum 
values), which corresponds to the long axis being along the 
line-of-sight. This result is the refore physically in tuitive, and 
consistent with the results in iRozo et al.l (|2007l ) for single- 
component isothermal ellipsoid lens models. The total vari- 
ation of otot with viewing angle for these typical triaxial 
models is a factor of ~ 2 at the lower mass scale, or as much 
as 4 in the higher mass case (where the flatter DM compo- 
nent contributes more). The 04/02 ratio reaches its maxi- 
mum value when a /a = 1 (maximum value) and b' /a — c/a 
(minimum value), which corresponds to the most flattened 
projected mass distribution. In the mixed case, the orienta- 
tion of the oblate stellar component is such that, regardless 
of (a'/i)dm, the top of the plot corresponds to a rounder 
projected stellar density because it is seen face-on, while the 
bottom of the plot corresponds to a flatter projected stel- 
lar density because it is seen edge-on, so that both otot and 
04/02 increase from top to bottom. As a result, the trends 
in both otot and 04/02 are not significantly modified from 
the triaxial case, except for a tilt in the contours which can 
be accounted for by considering the orientation of the oblate 
stellar component. 

We conclude that strong lenses drawn from a par- 
ent population with randomly-oriented non-spherical mat- 
ter distributions have a projected axis ratio distribution 
that is not consistent with that of a random sample of 
the population, because there are preferred orientations for 
the lensing systems. The degree of preference for particu- 
lar orientations can be factors of several depending on the 
intrinsic axis ratios, the relative contribution of the flatter 
DM component, and the number of images in the lenses 
being considered. While it is difficult to compare quanti- 
tatively with the conclusions of prev i ous works on "incli- 



nation bias" (e.g ., Mailer et al.lll997l ; iKeeton fc Kochanekl 
ll998l ; lRozoetal .1120071 ) due to the very different mass mod- 
els used, qualitatively we agree with the results that for an 
intrinsically non-spherical model, orientation effects are im- 
portant at this level. In particular, the trends we find re- 
garding otot and 04/02 fo r triaxial models agree well with 
those of IRozo et ail (|2007l ). 
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The preceding discussion has used results for fixed dark 
matter inner slope 7d m = 1- Generally, we find that the 
orientation bias trends in 04 are somewhat weaker for shal- 
lower 7dm = 0.5 models, and stronger for steeper 7dm = 1.5 
models, at the level of 10 to 30 per cent; whereas for 02 the 
degree of orientation bias is roughly the same for y^ m = 0.5 
and for 7<jm = 1.5. This may occur because using a shallower 
7dm makes the rounder stellar component more significant 
for the lensing properties, decreasing the overall flattening 
and hence decreasing the four-image cross-section. 

3.2 Shape 

Next we study orientation-averaged lensing cross-sections to 
see whether there are significant biases related to the intrin- 
sic shape of the galaxy. Unless explicitly noted, we again 
use our fiducial cusped models with 7d m = 1. Table [5] lists 
the total cross-sections and quad/double ratios of all seven 
shape models, for each magnification bias mode and for both 
mass scales. (Note that in some cases, the statistical errors 
are not identically zero but round to zero at the precision 
quoted in the table.) 

3.2.1 Total cross-section 

We first consider how the orientation-averaged unbiased to- 
tal cross-section varies with model shape for the lower mass 
model in Tabled When progressing from spherical to mod- 
erately and very oblate shapes, (crtot) decreases by 8 and 12 
per cent, respectively. When progressing from spherical to 
moderately and very prolate shapes, (crtot) increases by 6 
and 8 per cent, respectively. The triaxial model has (crtot) 
that is 3 per cent lower than in the spherical case. Finally, 
the cross-section for the mixed shape model is the same as 
that for the moderately oblate model, which suggests that 
the oblate stellar component is dominant in determining the 
lensing properties, even in the mixed model. Table [2] shows 
that these trends are only slightly enhanced for the higher 
mass model. We conclude that the orientation-averaged un- 
biased cross-section does not depend strongly on the intrin- 
sic model shape, even for significant changes in shape. In 
other words, while there is an orientation bias such that non- 
spherical lens galaxies are more likely to be aligned along 
the line-of-sight, there is no shape bias that would make 
non-spherical galaxies significantly over- or under-abundant 
in strong lens samples. 

Once we include magnification bias, the trends with 
model shape tend to be somewhat stronger, particularly 
when the survey flux limit is in the bright (steep) part of 
the luminosity function. However, even in the worst case, 
the variations in (crtot) with model shape are a maximum of 
~ 30 per cent, rather than ~ 10 per cent in the unbiased 
case. This is clearly not a selection bias that will drastically 
alter the distribution of model shapes for strong lensing sys- 
tems from the distribution of model shapes for all galaxies at 
that mass. Furthermore, while we have only tested specific 
sets of intrinsic axis ratios for the non-spherical shapes, the 
similarity of the cross-sections to the spherical case suggests 
that at least for axis ratios b/a and c/a equal to or greater 
than the cases we have tested, our conclusions still apply. 

With Table [2] we can consider how the shape of the 



galaxy affects the mass bias (at least for 7dm = 1; recall 
that Section [2~4l shows how the mass bias depends on 7d m for 
spherical models). For all model shapes, the unbiased total 
cross-section increases by a factor of ~ 3 from the lower to 
the higher mass scale. The difference is even larger when 
considering biased cross-sections, up to a factor of ~ 9 for 
cases with a flux limit in the steepest part of the luminosity 
function. We infer that mass bias is not very sensitive to 
galaxy shape (which is not surprising given that there is 
little dependence of (crtot) on shape). However, mass bias is 
sensitive to magnification bias such that in surveys with flux 
limits on the steeper part of the quasar luminosity function, 
the observed strong lensing systems will be preferentially 
skewed towards higher masses even more so than for those 
surveys with fainter flux limits. 

Our results com p are w ell with those in the literature. 
iKeeton fc Kochanekl (|l998T ) noted that there is a signifi- 
cant orientation bias for strong lensing by spiral galaxies, 
but found that the orientation-averaged cross-sections of 
disk+halo mo dels differ from thos e of spherical models by 
< 10 per cent. lHuterer et all |2005 f| used single-component 
projected isothermal ellipsoid models and found that, when 
the Einstein radius is fixed, the ellipticity affects lensing 
cross-sections at the few pe r cent level. One insta nce of 
a shape bias was found by lOguri fc Keetonl {2004), who 
found that (cr to t) increases with triaxiality for massive, gen- 
eralised NFW halos. However, they noted that the increase 
was driven by the cross-section for three-image lenses — an 
effect that was enhanced by the use of a single-component 
generalised NFW dark matter halo without a stellar com- 
ponent, such that the inner density profile was relatively 
shallow and hence quite prone to create three-image lens 
systems. Such models may or may not be reasonable for the 
massive halos of galaxy clusters, but are certainly not real- 
istic model_s_jOT_g_alMies ; Even so, it is interesting to note 
that lOguri fc Keetonl (|2004h found the net cross-section for 
two- and four- image lenses to be roughly independent of halo 
shape. Altogether, we conclude that there is general consen- 
sus that shape bias is not a strong effect in strong lensing 
by galaxies. 

Finally, we note that our conclusions about shape bias 
are not a strong function of mass model parameters such as 
7dm- This conclusion is expected given our agreement with 
other papers that use quite different mass models. However, 
as in Section \2. 41 and here, our conclusion regarding the de- 
gree of mass-based selection bias does depend on 7d m and 
the magnification bias mode. 



3.2.2 Quad/ double ratio 

Next, we consider the ratio of the orientation-averaged four- 
image and two-image cross-sections to understand the rel- 
ative abundances of quad and double lenses. Table [2] lists 



3 iKochanekl j 19961), IKeeton et al. (1997), Rusin &i Tegmarkl 
ll200ll) , and lChad (120031) " also computed lensing cross-sections for 
isothermal ellipsoid lens models with different ellipticities, but 
those studies did not isolate and discuss t he effects of shape on 
the total cross-section in the same way that lHuterer et all (2005) 
did. 
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Table 2. Lensing cross-sections averaged over all orientations, for each model shape and magnification bias mode, for the fiducial dark 
matter concentration and inner slope 7dm — It an d for both mass scales. 



Lower mass Higher mass 



Shape 


Bias mode 


(ftot) (arcsec 2 ) 


^4>/<<X2> 


(ctot) (arcsec 2 ) 


(ct4>/(o- 2 ) 


Spherical 


Unbiased 
Faintest, total 
Faintest, second 
Medium, total 
Medium, second 
Brightest, total 
Brightest, second 


0.53 ± 0.00 
1.08 ± 0.00 
0.72 ± 0.00 
1.97 ±0.01 
1.02 ± 0.01 
25.12 ± 0.22 
5.78 ± 0.10 


0.000 ± 0.000 
0.000 ± 0.000 
0.000 ± 0.000 
0.000 ± 0.000 
0.000 ± 0.000 
0.000 ± 0.000 
0.000 ± 0.000 


1.63 ±0.00 
4.18 ±0.01 
2.97 ±0.01 
8.34 ±0.02 
5.09 ±0.02 
163.52 ±0.85 
51.56 ±0.56 


0.000 ± 0.000 
0.000 ± 0.000 
0.000 ± 0.000 
0.000 ±0.000 
0.000 ± 0.000 
0.000 ± 0.000 
0.000 ±0.000 


Moderately oblate 


Unbiased 
Faintest, total 
Faintest, second 
Medium, total 
Medium, second 
Brightest, total 
Brightest, second 


0.49 ± 0.00 
0.98 ± 0.00 
0.64 ±0.00 
1.77 ±0.01 
0.87 ±0.00 
22.09 ± 0.09 
4.01 ± 0.02 


0.054 ± 0.001 
0.089 ± 0.001 
0.093 ± 0.001 
0.107 ±0.002 
0.132 ± 0.002 
0.251 ± 0.004 
0.611 ± 0.009 


1.46 ±0.01 
3.65 ±0.02 

2.47 ±0.01 
7.21 ±0.04 
4.08 ±0.02 

138.28 ±0.78 
36.16 ±0.23 


0.264 ± 0.004 
0.366 ± 0.005 
0.362 ± 0.005 
0.406 ± 0.006 
0.429 ± 0.006 
0.645 ±0.009 
1.055 ±0.013 


Very oblate 


Unbiased 
Faintest, total 
Faintest, second 
Medium, total 
Medium, second 
Brightest, total 
Brightest, second 


0.47 ±0.00 
0.93 ± 0.00 
0.60 ±0.00 
1.66 ± 0.01 
0.81 ± 0.00 
20.64 ±0.10 
3.90 ± 0.03 


0.110 ±0.002 
0.168 ± 0.002 
0.173 ± 0.002 
0.198 ± 0.003 
0.239 ± 0.003 
0.448 ± 0.006 
0.957 ±0.013 


1.48 ±0.01 
3.73 ±0.03 
2.51 ±0.02 
7.38 ±0.06 
4.15 ±0.03 
143.55 ± 1.14 
37.46 ±0.33 


0.540 ±0.010 
0.684 ±0.012 
0.662 ±0.011 
0.742 ±0.013 
0.750 ±0.012 
1.077 ±0.017 
1.465 ±0.021 


Moderately prolate 


Unbiased 
Faintest, total 
Faintest, second 
Medium, total 
Medium, second 
Brightest, total 
Brightest, second 


0.56 ± 0.00 
1.14 ±0.00 
0.74 ±0.00 
2.06 ± 0.00 
1.01 ± 0.00 
25.61 ± 0.09 
4.78 ± 0.03 


0.010 ± 0.000 
0.020 ± 0.000 
0.021 ± 0.000 
0.025 ± 0.001 
0.032 ± 0.001 
0.064 ± 0.001 
0.166 ± 0.003 


1.77 ±0.00 
4.49 ±0.01 
3.07 ±0.01 
8.92 ±0.02 
3.07 ±0.01 
172.87 ±0.56 
45.60 ±0.23 


0.035 ±0.001 
0.057 ±0.001 
0.058 ±0.001 
0.065 ±0.001 
0.058 ±0.001 
0.109 ±0.002 
0.203 ±0.003 


Very prolate 


Unbiased 
Faintest, total 
Faintest, second 
Medium, total 
Medium, second 
Brightest, total 
Brightest, second 


0.57 ± 0.00 
1.15 ± 0.00 
0.74 ±0.00 
2.07 ±0.01 
1.00 ± 0.00 
25.51 ± 0.10 
4.59 ± 0.03 


0.022 ± 0.000 
0.040 ± 0.001 
0.042 ± 0.001 
0.049 ± 0.001 
0.062 ± 0.001 
0.119 ± 0.002 
0.304 ± 0.006 


1.83 ±0.01 
4.60 ±0.01 
3.13 ±0.01 
9.11 ±0.03 
5.19 ±0.02 
174.92 ± 0.64 
45.29 ±0.25 


0.078 ±0.001 
0.117 ±0.002 
0.117 ±0.002 
0.131 ±0.002 
0.142 ±0.002 
0.210 ±0.003 
0.370 ± 0.005 


Triaxial 


Unbiased 
Faintest, total 
Faintest, second 
Medium, total 
Medium, second 
Brightest, total 
Brightest, second 


0.51 ± 0.00 
1.03 ± 0.00 
0.66 ± 0.00 
1.86 ± 0.00 
0.90 ±0.00 
23.07 ±0.06 
3.96 ± 0.02 


0.042 ± 0.000 
0.073 ± 0.001 
0.077 ± 0.001 
0.089 ± 0.001 
0.113 ± 0.001 
0.218 ± 0.002 
0.636 ± 0.006 


1.55 ±0.00 
3.86 ±0.01 
2.60 ± 0.01 
7.62 ±0.02 
4.27 ±0.01 
144.19 ±0.44 
35.99 ±0.14 


0.176 ±0.002 
0.260 ±0.002 
0.261 ± 0.002 
0.292 ± 0.002 
0.321 ±0.003 
0.490 ± 0.004 
0.983 ±0.007 


Mixed 


Unbiased 
Faintest, total 
Faintest, second 
Medium, total 
Medium, second 
Brightest, total 
Brightest, second 


0.49 ± 0.00 
0.98 ± 0.00 
0.63 ± 0.00 
1.77 ±0.00 
0.86 ± 0.00 
22.24 ±0.06 
3.94 ±0.02 


0.059 ± 0.001 
0.099 ± 0.001 
0.103 ± 0.001 
0.118 ± 0.001 
0.148 ± 0.001 
0.281 ± 0.003 
0.752 ± 0.007 


1.47 ±0.00 
3.69 ±0.01 
2.49 ±0.01 
7.28 ±0.02 
4.10 ±0.01 
139.49 ±0.46 
35.64 ±0.14 


0.220 ±0.002 
0.316 ±0.003 
0.317 ±0.003 
0.354 ±0.003 
0.385 ± 0.003 
0.582 ±0.005 
1.103 ±0.008 



(a4}/{<72) for all shape models, for both mass scales. Natu- 
rally, in the spherical case this ratio is identically zero. For 
oblate models, the abundance of four-image lens systems 
can be quite significant, particularly when the flux limit is 
in the steep part of the luminosity function: for example, for 



the very oblate, higher mass model the quad/double ratio 
ranges from ~ 0.5 in the unbiased case to ~ 1.5 for some 
of the magnification bias modes. Prolate systems, in con- 
trast, have fairly low quad/double ratios, reaching a maxi- 
mum of 37 per cent. The trend with model shape is sensi- 
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ble: for the moderately and very oblate cases the smallest 
axis ratio (c/o)dm is (respectively) 0.50 and 0.36, whereas in 
the moderately and very prolate case the smallest axis ratio 
(fo/a)dm = (c/a)dm is (respectively) 0.71 and 0.60. Thus, the 
maximum projected ellipticity is lower in the prolate case, 
so there are fewer four-image lenses. 

The trend with magnification bias is due to the fact that 
four-image systems have preferentially higher magnifications 
than two-image systems. Moreover, four-image lenses often 
have two or more particularly bright images, which is why 
the quad/double ratio is especially high if we compute mag- 
nification bias using the second-brightest image (see also 
Fig. [1]). In the triaxial case, the quad/double ratio is par- 
ticularly dependent on how magnification bias is treated, 
ranging from 0.04 to 1.00. 

Generically, for a given shape model we find that 
(er4)/(<T2) is larger for our higher mass model, because at 
the higher mass scale the dark matter component (which 
is flatter than the stellar component) plays a more impor- 
tant role. Our conc lusion is seemingly at odds with that of 
IMoller et all l|2007h . who found that the four-image systems 
in their simulations are preferentially at lower velocity dis- 
persio n than the two-image systems. However, IMoller et al.l 
l)2007h also found that the four-image systems have lens 
galaxies with a more pronounced disk component, in which 
case the appropriate comparison is not at fixed DM halo 
shape. Thus, there is no real disparity between our re- 
sults and theirs, since we find that oblateness (as would oc- 
cur with a pronounced disk component) strongly increases 

(o-4>/(cr 2 ). 

Our results are compatible with the conventional wis- 
dom that the quad/double ratio in strong lensing offers a 
valuable probe of the density shapes in the inner r egions 
of galaxies (e.g.. lRusin fc Tegmarkll200ll ; IOgurill2007l ). But 
magnification bias also plays a very im portant role (s ee also 
iHuterer et"alll2005l : IMoller et alJl2007ft . In particular. lOeuril 
<|2007f ) suggests that magnification bias may explain why 
there is a higher fraction of four-image lenses in radio sur- 
veys than in optical surveys. He suggests that the source 
luminosity function is steeper in radio surveys than in op- 
tical surveys (at least, in optical surveys that probe deeper 
than the break in the quasar luminosity function), which 
enhances magnification bias and (as we have seen) increases 
the quad/double ratio in radio surveys. 

3.2.3 Three-image cross-section 

The final shape-related question we address is which shapes 
give rise to three-image lens systems. They can only occur 
in cases of such high projected ellipticity that the tangential 
caustic protrudes beyond the radial caustic. For the lower- 
mass model, we find three-image systems only for the very 
oblate model, with a cross-section that depends sensitively 
on the dark matter inner slope. This dependence may arise 
from the fact that our dark matter is flatter than the stel- 
lar components, and models with steeper dark matter inner 
slope are more dark matter-dominated in the region where 
strong-lensing is important. Thus, for example, the three- 
image cross-section (unbiased) for this very oblate lower 
mass model averages to zero for 7a m = 0.5 and 7dm = 1, 
and is small but non-zero for 7dm = 1-5 (but (0-3) is still 
just 0.3 per cent of (otot))- Magnification bias can enhance 



the relative abundance of three-image systems by up to a 
factor of ~ 10 depending on the survey flux limit (with the 
greatest enhancement when the flux limit is on the bright 
end). Naturally these systems can only occur in the oblate 
case when the galaxy is seen close to edge-on, which is a 
fairly extreme orientation-related selection bias inherent to 
such systems. 

For the higher-mass model, three-image lens systems 
occur for both the moderately and very oblate cases (pre- 
sumably because the higher-mass model is inherently more 
dark matter-dominated and therefore flatter); and there is 
again a strong trend with 7d m - For the very oblate shape, the 
unbiased three-image lens system abundance ((<73)/(<7tot)) 
ranges from 0.7 to 8 per cent as 7dm changes from 0.5 to 1.5, 
a factor of 10 increase. For the most extreme magnification 
bias mode (4L* with the total magnification) , the increase is 
from ~ 3 per cent at 7dm = 0.5 to ~ 20 per cent at 7dm = 1-5 
(a factor of 7). An observational scarcity of three- image lens 
systems with group-mass lenses can therefore constrain the 
abundance of strongly oblate matter distributions (but see 
below). For this mass scale, the triaxial and mixed models 
also give rise to three-image lenses, but with a relative abun- 
dance that is below 1 per cent for all 7dm and magnification 
bias modes. 

We contra s t th is result with the findings in 
lOguri fc KeetorJ (|2004l ). who use single-component gener- 
alised NFW profiles with an axis ratio distribution (rather 
than fixed shapes as used here) to study the trends of 
(o"3)/(o"tot) with 7dm- That work suggests that the three- 
image abundance increases as 7d m gets shallower, which 
seems to contradict our trend. The difference between 
these two results stems from a difference in modeling as- 
sumptions: in their case, without a fixed stellar component 
creating a lower bound to the central density, decreasing 
7d m strongly affects the (lens-plane) inner critical curve 
(which is very sensitive to the central slope of the den- 
sity profile), shrinking the (source-plane) radial caustic. 
However, the tangential caustic does not shrink as much, 
because it depends primarily on the projected ellipticity. 
As a result, decreasing 7d m actually increases the area that 
is inside the tangential caustic and outside the radial one, 
where three-image systems originate. In contrast, when 
we decrease 7d m the radial caustic does not shrink very 
much, since we have a fixed stellar component. However, 
the stellar component becomes much more dominant, and 
since we have modeled it as being rounder, this means the 
tangential caustic will shrink. The result is that the area 
inside the tangential caustic and outside the radial one, 
and therefore the three-image lensing system abundance, is 
reduced. 

This comparison suggests that the predicted abundance 
of three-image strong lenses is quite dependent on modeling 
assumptions, particularly on the ratio of stellar to dark mat- 
ter mass and the assumed shapes for the two components. 
Another potential concern in studying the statistics of three- 
image lenses is that a four-image lens could be misclassified 
as a three-image system if one of the images is too faint 
to be detected. Since typically (0-4) S> (03}, any small rate 
of contamination from misclassified four-image lenses could 
corrupt the three-image sample. In practise, the frequency 
with which such errors may occur will depend sensitively on 
the galaxy model and on the survey parameters, so a detailed 
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exploration is beyond the scope of this paper. For now, we 
conclude that the abundance of three-image lenses contains 
interesting information about the abundance of very oblate 
galaxies, but there may be practical challenges in extracting 
that information from observed lens samples. 



3.3 Magnification bias 

As an aside, we now investigate whether the conversion from 
unbiased to the various biased cross-sections can be quanti- 
fied with a single number for each magnification bias mode, 
independent of mass, shape, dark matter inner slope, or the 
number of images. This would allow considerable simplifica- 
tion, not only in the presentation of our results, but also 
in future calculations. We again refer to the orientation- 
averaged cross-sections in Tableful 

First, we find that for 7d m = 1 (as in the table), and 
for a given mass scale, the conversion from unbiased to bi- 
ased otot is relatively independent of shape. This statement 
is true within ~ 2 per cent for the fainter flux limits, or 
within < 10 per cent for the brightest flux limit. However, 
the conversion does depend strongly on 7d m and mass. In 
retrospect, this is relatively easy to understand, because the 
biased cross-section depends not just on the lensing caus- 
tics, but also on the lensing magnification as a function of 
position, which in turn depends on how the projected den- 
sity is distributed within the Einstein radius. Changing the 
dark matter inner slope changes both the caustics and the 
magnification distribution, and not necessarily in the same 
way. Likewise, the different mass scales have different bal- 
ances of stellar and dark matter components, each with its 
own density profile, so the conversion from unbiased to bi- 
ased cross-sections is not the same for our lower and higher 
mass scales. 

Furthermore, because the four-image cross-section de- 
pends on the projected shape (Section 13. 2p . the conversion 
of the quad/double ratio from the unbiased to the biased 
cases depends not just on 7d m and the mass, but also on the 
shape. We therefore conclude that for the mass models used 
here, there is no simple conversion from unbiased to biased 
cross-sections that would allow us to simplify the presenta- 
tion of our results. Selection biases in <T to t and 0-4/(72 with 
the form of the mass density profile may depend strongly on 
the survey detection strategy and flux limit. Consequently, 
we continue to present results for the different magnification 
bias modes throughout the rest of this paper. 
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Figure 3. Variation of total, orientation-averaged lensing cross- 
section with 7d m for the lower (left) and higher (right) mass scales 
for different shape models as labeled on the plot, with the mag- 
nification bias mode indicated as in Fig. [T] 
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4 INNER SLOPE 

We now examine the role of the galaxy density profile in 
more detail. We study how the dark matter inner slope 7d m 
affects orientation-averaged cross-sections (Section I4.1|l as 
well as lens image separation distributions (Section 14. 2[) . In 
Section 14.31 we extend the discussion of profile-dependent 
selection biases to consider the effect of adiabatic contrac- 
tion on the strong lensing cross-section (for spherical models 
only). 



Figure 4. Variation of the orientation-averaged quad/double ra- 
tio with 7d m for the lower (left) and higher (right) mass scales for 
different shape models as labeled on the plot, with magnification 
bias mode indicated as in Fig. [T] 

4.1 Inner slope 

In the sections that follow, we typically use orientation- 
averaged cross-sections. While in sections that dealt with 
orientation and halo shape, we have introduced the notation 
(o" t ot) to highlight the difference between cross-sections for a 
single orientation versus orientation-averaged cross-sections, 
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we now revert to simply using a to t for orientation- averaged 
cross-sections since they are now our sole focus. 

Fig. [3]shows the variation of the total cross-section with 
the dark matter inner slope 7dm for both mass scales, for four 
of the seven shapes and all seven magnification bias modes. 
In order to highlight the variation with 7d m (and to fit all 
results for a given shape in a single panel), we normalise 
<Ttot for each magnification bias mode by the corresponding 
values when 7d m = 1 (given in Tabled • Statistical error bars 
are shown on the plot, but they are generally smaller than 
the sizes of the points themselves due to the large number of 
Monte Carlo realisations. Fig. [4] shows analogous results for 
the quad/double ratios (once again, normalised by the values 
at 7d m = 1 given in Table [2}. We again show different shape 
models, although we eliminate the spherical model since its 
four-image cross-section is identically zero. 

The total lensing cross-section for the lower mass model 
(Fig. [3j left column) shows several interesting trends with 
7dm- First, as we expect, a to t increases with 7d m because 
the projected mass density is higher in the inner parts for 
steeper 7dm- Second, all the curves for a given mass scale are 
nearly the same, with er t ot decreasing by 20 to 30 per cent 
when going from 7d m = 1 to 7a m = 0.5, and increasing by a 
factor of ~ 3 when going from 7d m = 1 to 7a m = 1.5. The 
variation within the quoted ranges depends only slightly on 
the magnification bias mode, and is nearly independent of 
the shapes considered here. For the higher mass scale (Fig.|3l 
right column) , we likewise see that the results are consistent 
for different shapes. At this mass scale, <r t ot decreases by 50- 
70 per cent going from 7d m = 1 to 7d m = 0.5, and it increases 
by factors of 6-10 going from 7d m = 1 to 7d m = 1-5 (where 
the quoted range encompasses the different magnification 
bias modes). The reason the dark matter inner slope causes 
more significant selection bias for the higher mass model 
is that the flatter dark matter is more prominent in that 
model. 

To understand the strong variation with 7d m further, 
we explore a finer grid in 7d m in Section f5. li b elow , and we 
compare the variation in the cross-section with the variation 
in the projected dark matter fraction. However, it is appar- 
ent from paper I figure 6 that the projected mass density 
(and therefore the lensing properties) is dominated on small 
scales by the stellar component for 7d m = 0.5 and 7d m = 1, 
whereas the dark matter component has a significant contri- 
bution everywhere for 7dm = 1.5. Consequently, the cross- 
section variation is quite likely related to the change in the 
projected dark matter fraction. 

Next, we consider the quad/double ratio for the lower 
mass model in the left column of Fig. [4] For all shapes 
shown here, this ratio is roughly exponential in 7dm, 04/02 oc 
exp (0.857dm), once again roughly independent of model 
shape and with consistent trends for all magnification bias 
modes. Moving to the higher mass scale in the right column, 
we see rather different results: the quad/double ratio in- 
creases from 7dm = 0.5 to 7dm = 1, but then approximately 
flattens out from 7d m = 1 to 7dm = 1.5. There is actually a 
slight decrease in 4/ '02 for the unbiased cross-section, and a 
slight increase for the significant magnification bias modes. 
This result involves a combination of several effects, the sizes 
of which approximately cancel, and which are complicated 
by the fact that we are working with orientation-averaged 
quantities. First, examination of the mock lens systems sug- 



gests that the radial caustic moves strongly outwards for 
7dm = 1-5 for the higher mass scale (leading to the very large 
increase in to t shown in Fig. [3}. The tangential caustic in- 
creases in size somewhat, because it depends predominantly 
on the projected shape (which is flatter at high 7dm due to 
the dominant DM component), so while 04 increases, so does 
02 ~ 0"tot — 04, thus leading to a flattening out of 04/02. Be- 
cause the increase in to t is not nearly as much for the lower 
mass model when going from 7d m = 1 to 1.5 (due to the 
lower dark matter fraction), 04/02 does not level out as for 
the higher mass model. The trend with magnification bias is 
also understandable: for more significant magnification bias 
(higher flux limit and steeper luminosity function), the inner 
regions where 4-image systems dominate become more im- 
portant, which is why for the higher mass model the biased 
04/02 can still increase slightly when going from 7dm = 1 to 

1.5 even though the unbiased qua d ratio decreases. 

Several previous st u dies (e.g., iPorciani fc Madaull2000l; 



Keeton fc Madaul l200ll: IWvithe et all l200ll : iLi fc Ostrikerl 
20021 ; iHuterer fc Mall2004n ~ used spherical single-component 



mass models to examine how the lensing cross-section 
depends on the inne r slope of the density profile. 
iKeeton fc Madaul l|200ll ) found in particular that the cross- 
section variation could be explained mainly by changes in 
the fraction of the total mass that lies within some fiducial 
radius. The simple scaling with the "core mass fraction" is 
probably broken in our models by the presence of the stel- 
lar component, although we do examine whether there is 
some analogous scaling in Section 15.11 below. There is gen- 
eral consensus that the total cross-section can increase fairly 
dramatically as the inner density profile steepens, although 
we believe our two-component models better isolate the ef- 
fects of the dark matter slope alone. 

In conclusion, there is a strong selection bias in the total 
cross-section with 7d m that can amount to factors of a few 
(particularly within the upper end of the allowed range); 
and this bias is remarkably insensitive to the model shape 
and only mildly sensitive to the magnification bias (survey 
flux limit and detection strategy). The selection bias with 
7dm in the quad/double ratio can amount to several tens of 
per cent, and is likewise relatively insensitive to the model 
shape and magnification bias mode. 

In this section, we have focused on quantities averaged 
over orientation. In Appendix [A] we indicate how the vari- 
ation of cross-section with 7d m depends only mildly on ori- 
entation. This result is relevant only when a special subset 
of strong lenses are selected based on their projected shape. 



4.2 Image separation distributions 

We next consider the effects of the dark matter inner slope 
7dm on the image separation distribution, and consider 
whether this is a magnification bias-dependent quantity. 
This calculation is done separately in the spherical and the 
non-spherical cases. 



4-2.1 Spherical models 

We begin with cusped, spherical density profile models, 
which lead only to strong lensing systems with two images. 
We first consider what to expect for the image separation 
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Figure 5. Mean (points) and width (errorbars) of the image 
separation distribution for the spherical shape model, in units 
of 2i?, c ; n , as a function of 7d m - Points are horizontally offset 
for clarity, but all measurements are with the values of 7d m = 
{0.5, 1, 1.5}. Line colours and types indicate the magnification 
bias mode, as described in the caption for Fig. [T] 



distribution based on general considerations. By definition, 
the Einstein radius 7? c i n for a circularly-symmetric projected 
mass distribution is the radius at which the lensed image is a 
ring because the source is located precisely behind the lens. 
In that case, the image separation is AO — 2R c i n . Conse- 
quently, we expect that regardless of the form of the density 
profile, it must approach 2_R e i n , because the sources near the 
source plane origin produce precisely that image separation. 

The simplest comparison we can make is against the 
very simple image separation distribution for the SIS density 
profile, which is a delta function at 2J? c i n : 



psis(A0) = 8(A0 - 2.Rei n ). 



(•>) 



We ask, then, how the complexity of our mass models causes 
the image separation distribution to deviate from this simple 
form. 

Fig. [5]shows the mean image separation in units of twice 
the Einstein radius as a function of 7dm, for all magnification 
bias modes. The standard deviation of the image separation 
distribution is shown as an errorbar. There are a few conclu- 
sions we can draw from this plot. First, the image separation 
distribution depends on 7d m almost exclusively through its 
dependence on 7? e in- Second, despite the complexity of the 
mass models, to lowest order the distributions are remark- 
ably close to the distributions for a pure SIS, very narrow 
and with an average that is only ~ 2-3 per cent below 2R c i n . 
Finally, there is a minor (for practical purposes) but nonzero 
dependence on the magnification bias mode. The mean of 
the image separation distribution increases with more sig- 
nificant magnification bias, approaching 2R e i u from below, 
because the magnification is highest for a source near the 
origin, where A9 is precisely 2i? c i n . 



We have quantified the image separation statistics using 
the usual mean and standard deviation. However, the dis- 
tributions p(A0) are not remotely Gaussian or even peaked 
at (AO); rather, they are steeply increasing functions of A6 
that are truncated abruptly at A6 = 2i? e i n . We use the de- 
flection angle profile a(R') in paper I figure 7 to explain the 
shape of this distribution, and why the mean is so close to 
the SIS result. As shown there, for the three cusped models, 
a{6) is nearly flat and equal to 7? e in from R' ~ i? c in/4 until 
well above i? e in- Consequently, the image separations should 
be below 2J? e in for those small subset of images that appear 
to be near the origin in the lens plane, but approach 2i? c i n 

for the majority of the images. 

Our results are consistent with those of lOguri et al.l 

(|2002l ). who study image separations for spherical, single- 
component, generalised NFW models, and also find that it 
is close to, but slightly below, 2R c i n . In light of that result, 
ours is not surprising, since the superposition of a stellar 
component with a generalised NFW model tends to yield a 
composite model that is closer to isothermal near i? e i n and 
therefore more likely to give quasi-isothermal image separa- 
tion distributions. 

Thus, for spherical models, the changes in the image 
separation distribution due to changes in density profile pa- 
rameters (mass, 7dm) arise almost exclusively through the 
change in the Einstein radius. This is a remarkably simple 
result given the complexity of the mass models used. We 
note, however, that while the physical scaling of the image 
separation distribution is simple, there may be observational 
selection biases due to resolution limits (more likely to miss 
lenses with smaller R c i n ) or due to finite apertures for the 
image search (more likely to miss lenses with larger R c in)- To 
estimate the effects of selection bias in any given survey, it 
will be necessary to supplement the physical selection effects 
we study here with an analysis of observational selection ef- 
fects. 



4-2.2 Non-spherical models 

For non-spherical models, we have defined the image separa- 
tion as the maximum separation between any pair of images 
in a given image system, regardless of the total number of 
images. We have then calculated the mean and width of the 
image separation distribution as a function of 7dm for both 
mass scales, with different extreme configurations in the pro- 
jected shapes. These results are shown for the unbiased case 
in Fig. [6] normalised by twice the spherical Einstein radius 
(i.e., we use the same normalisation factors that went into 
Fig. [5}. As for the spherical case, the effects of magnification 
bias (not shown) are subtle. 

There are a few conclusions we can draw from this plot. 
First, as expected, the mean of the image separation distri- 
bution scales with 7d m in nearly the same way as the Ein- 
stein radius in the spherical case, so the curves are nearly 
flat. Second, the difference between the means and widths 
of p(A8) for various shape configurations can be understood 
simply in terms of the projected shape. For example, for the 
very oblate and very prolate shapes, the solid lines (edge- 
on and end-on, respectively) give the highest mean A6 be- 
cause they give a larger projected mass density for a given 
(azimuthally-averaged) separation from the centre of mass 
than in configurations represented with dashed lines (face- 
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Figure 6. Mean (points) and width (errorbars) of the unbiased 
image separation distribution, relative to 2i? c ; n for the spheri- 
cal shape model, as a function of 7d m - The left and right panels 
show the lower and higher mass scales, respectively. Points are 
horizontally offset for clarity, but all measurements are with the 
values of 7(j m = {0.5, 1, 1.5}. The top panels show the very oblate 
model shape, where the solid line = edge-on, dotted line = in- 
termediate, and dashed line = face-on configuration. The middle 
panels show the very prolate model shape, where the solid line 
= end-on, dotted line = intermediate, and dashed line = side-on 
configuration. The bottom panels show the triaxial model shape, 
where the solid line = viewed with line-of-sight along a, dotted 
line = viewed along b, and dashed line = viewed along c. 



on and side-on, respectively), consistent with the larger to- 
tal lensing cross-section. Because the projected ellipticity 
is higher for the solid line configuration in the oblate case 
(edge-on) and for the dashed line in the prolate case (side- 
on), the width of the image separation distributions in those 
cases is larger (depending on where the source is located in 
the region where lensing occurs, very different image con- 
figurations can occur). In the triaxial case, these arguments 
imply that the mean image separation is highest when view- 
ing along the major axis (solid lines), but the width is high- 
est when viewing along the intermediate axis (dotted lines) . 
When considering the image separation distribution derived 
from many lenses drawn from galaxies of a given shape with 
random orientations, we must average the image separa- 
tion distribution for each orientation weighted by the cross- 
section CTtot for that orientation. We have confirmed that 
the differences between (AO) for different orientations rela- 
tive to the spherical Einstein radius can be accounted for by 
computing ii e i n from the monopole deflection (equation 26 

Of Paper I) as ao(-Rcin) = -Rein- 

Our finding that non-spherical models have significantly 
broader image separation distributions (~ 10 per cent rather 
than a few per cent), but mean values that are not strongly 
different from the sph e rical c ase, is consistent with the re- 
sults in iHuterer et al.1 |2005) , who tested this effect using 
single-component singular isothermal ellipsoid models. The 



fact that the image separation distribution is quite close to 
that for a spherical isothermal model even for these highly 
non-spherical and non-isothermal models is quite remark- 
able. The implication is that there may be an observational 
selection bias that is sensitive to 7dm, depending on the an- 
gular resolution and aperture limits for a particular survey. 
For example, an inability to find large-separation lens sys- 
tems would tend to counteract the physical selection bias 
towards galaxies with steeper dark matter inner slope (7d m ), 
particularly at low redshift. 

4.3 Adiabatic contraction 

As discussed in Paper I, the need for adiabatic contraction 
of dark matter halos in res ponse to the adiabatic infall of 
baryons to form a galaxy dYounelll980l; iBlumenthal et al.1 
1 19861 ; iGnedin et af1l2004l ; ISellwood fc McGaughll2005r i is un- 
clear for early- type galaxies. In this subsection, we test 
the effects of adiabatic cont raction explicitly, fo r sphe rical, 
cusped 7dm = 1 models onlv. lKochanek fc White! (|200ll ) pre- 
viously examined the effect of AC on the s trong l ensing im- 
age separation distribution, and iKeetonl (|2001al ) included 
AC when computing lens statistics, but neither study con- 
sidered it as a selection bias as we are. 

The effects of AC on the projected density profiles are 
shown in paper I figures 4-7; it increases the dark matter 
density on scales below ~ 10 kpc. Due to the increased 
density in the inner regions, we expect an increase in the 
strong-lensing cross-section. Our purpose is to determine the 
significance of this increase compared to the increase when 
changing 7 dm = 1 to 7 dm = 1.5 without AC. 

For the lower mass model, the total lensing cross-section 
Otot increases by a factor of 2.6 when going from 7d m = 1 to 
7dm = 1.5, or slightly less than that for the higher magni- 
fication bias modes. In contrast, when going from 7d m = 1 
without AC to includi ng AC, the cross-section increases by a 
factor of ~ 2.6 for the lBlumenthal et al.l (Il986h prescription, 
or a factor ~ 2.0 for the Gnedin et al.1 (|2004l ) prescription. 
These factors are reduced to ~ 2.0 and ~ 1.5 for the higher 
magnification bias modes. 

For the higher mass model, the unbiased otot increases 
by a factor of 9 when going from 7d m = 1 to 7d m = 1.5, 
or as little as 6 for the higher magnification bias modes. 
In contrast, when going from 7d m = 1 without AC to in- 
clu ding AC, otot incr e ases b y a factor ~ 3.5 for the origi- 
nal HhHnenthaTeF]^ §wM) prescription, or a factor ~ 2.3 
for the Gnedin et al.1 (|2004| ) prescription. These factors are 
again reduced for the higher magnification bias modes. For 
comparison, previous studies of the effect of baryonic cool- 
ing on strong lensing by clusters also foun d that the cross- 
section increases by a factor of a few (e.g.. | Puchwein et al.1 
120051 ; iRozo et altood IWambsganss et alJl20ojT ~ 

We conclude that AC has typically a less significant 
effect on the strong lensing cross-section than increasing 
7dm from 1 to 1.5 for both mass scales. Nonetheless, the 
increases in the strong lensing cross-section due to AC are 
non-negligible. There are several implications of this finding. 
First, if AC happens in typically all galaxies of these types 
(i.e., the assumptions behind it are valid and the changes in 
profile are preserved through mergers), it suggests that the 
numbers of strong lenses in upcoming surveys will be signif- 
icantly higher than predicted given a simple two-component 
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model that does not include AC. Second, if AC only oc- 
curs in massive, early-type galaxies with particular shapes 
or formation histories, then the ones for which AC does oc- 
cur are more likely to be strong-lensing systems. Thus, for a 
given mass system, there may be a selection bias favouring 
galaxies or groups depending on their formation history. 

Improved simulations, better understanding of galaxy 
formation, and observational results that confirm, contra- 
dict, or refine these AC models will be useful for understand- 
ing selection biases due to AC. A more detailed exploration 
is beyond the scope of this paper. 



5 CONCENTRATION 

We now allow more model freedom than in the previous 
section, varying both the dark matter inner slope 7d m and 
the dark matter halo concentration Cdm on a grid (spherical 
model only) in Section 15.11 We also consider a more lim- 
ited set of changes in dark matter halo concentration in the 
triaxial case in Section \5. 2 1 

5.1 Concentration versus inner slope 

Here we explore the trade-off between steepening the dark 
matter inner slope versus increasing the concentration (with 
fixed stellar component), both of which are expected to in- 
crease the lensing cross-section due to the enhanced density 
of dark matter in the very inner parts of the galaxy. 

First, we explore this degeneracy with a relatively fine 
grid in dark matter inner slope jdm and concentration Cdm, 
for spherical cusped models only. Our grid has 10 values 
of 7dm and 8 values of Cd m for each mass scale. In accor- 
dance with our finding in Section 14.11 that the increase in 
total cross-section with 7d m for slopes steeper than unity is 
extremely rapid at our fiducial concentration, we choose 10 
values of 7d m from to 2, distributed evenly in exp(7d m ). 
To define the grid in concentration, we use the result from 
JV-body simulations that the distribution of Cdm values is 
roughly lognormal around the fiducia l concentration, with 
scatter of 0.15 dex l|Bullock et alj|20 01). Our 8 values of Cdm 
range from — 3<j to +4a, i.e., from Cd m = 3.0 to Cd m = 33.4 
for the lower mass scale, and from Cdm = 2.0 to Cdm = 22.3 
for the higher mass scale, in steps of la. (The grid is skewed 
towards higher values because of our expectation that the 
cross-section will get quite large in those trend we 

want to track). When we change 7dm and/or Cdm, we nor- 
malise the profile by fixing Miso and allowing po to vary, as 
discussed in Paper I. 

For some of the grid points at high 7d m and/or Cdm, 
the Einstein radius is significantly larger than for our test 
case (with 7d m = 1.5 and fiducial concentration) that we 
used to determine adequacy of the box size in Paper I. To 
avoid box size effects for those grid points we increase the 
box size by a factor of 3/2 (lower mass scale) or 2 (higher 
mass scale), while maintaining the pixel size. We expect any 
remaining box size effects to reduce the cross-section for the 
steepest slopes and highest concentrations, so any trends we 
see with these parameters could be slightly underestimated. 
For grid points at lower 7d m and/or Cd m than our test case, 
the Einstein radius is smaller than the 7dm = 0.5, fiducial 
concentration case for which we tested resolution. However, 



we find that it is at most 15 per cent smaller, because of the 
unchanging stellar component which means the combined 
profile cannot get arbitrarily small, and since our resolution 
choice was conservative in the spherical case, we do not in- 
crease the resolution any further. 

For each of the 80 grid points, we compute the unbiased 
and biased lensing cross-sections; since the models are spher- 
ical, there is only one cross-section to consider (<T to t = 02, 
(73 = o"4 = 0). Fig. [7] shows contour plots of the base-10 
logarithm of the lensing cross-section for both mass scales. 
For illustration purposes, we show the unbiased cross-section 
only, and describe the effects on the biased cross-sections in 
the text. For the lower mass scale (left panel), the total 
variation of the unbiased cross-section is a factor of 200; for 
the higher mass scale (right panel), a factor of 3 200. This 
difference in dynamic range for the two mass scales is be- 
cause the higher mass scale was already more dark matter- 
dominated than the lower mass scale even for 7d m = 1 and 
the fiducial concentration, so increasing either one or both of 
those values has a more dramatic effect on the cross-section. 
Much of this strong variation with cross-section occurs in a 
fairly small range of 7d m ^ 1-3, and the very highest cross- 
sections are only achieved for concentrations that are several 
a above the median. Nonetheless, the enormous variation in 
cross-section across the parameter space suggests that the 
observed distribution of concentration and inner slopes from 
strong lensing systems, psi.(7dm, Cdm), may be quite biased 
compared with the intrinsic distribution (at fixed mass). In- 
deed, in Section \S . 1 1 below we construct an example to illus- 
trate the bias in the distribution of 7d m and Cd m that may be 
observed in strong lensing surveys. We find that the distribu- 
tion can be both shifted and broadened by the lensing selec- 
tion bias. These effects produce a significant (factor of sev- 
eral) enhancement in the abundance of high-concentration 
halos that, according to the intrinsic distribution, should 
constitute only ~ 1 per cent of the population. 

Next, we look at slices through the grid. For the fiducial 
concentration itself, the total variation in the cross-section 
with 7dm is a factor of 40 (lower mass scale) or nearly 600 
(higher mass scale) going from 7d m = to 7d m = 2. For the 
fiducial inner slope (7dm = 1), as the concentration varies 
from — 3<t to +4<r, the cross-section varies by a factor of 
five (lower mass) to several tens (higher mass). Thus, we 
can conclude that (a) the effects of varying concentration 
and inner slope of the dark matter are significant, and (b) 
those effects are (not surprisingly) stronger in systems where 
the dark matter is more dominant. In either case, it seems 
clear that efforts to derive the underlying distribution of 7d m 
and Cdm values using strong lensing-selected systems must 
account for this strong bias in the cross-sections with dark 
matter inner slope and concentration. 

We also find that the trends of the cross-section with 
7dm and Cdm do not depend too significantly on the mag- 
nification bias modes. The total variation across the grid is 
somewhat smaller for the biased cases than for the unbiased 
case, because for the steeper slope models the cross-section 
does not increase as sharply as in the unbiased case. 

We further explore the degeneracy between 7d m and 
Cdm by considering, for all models on the grid, the Einstein 
radius _R e in and the (projected) dark matter fraction / dm 
within i? e i n . We expect both to be connected to the lensing 
cross-section cr t ot in some non-trivial way. Fig. [8] shows the 
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Figure 7. Contour plots of the base-10 logarithm of the unbiased cross-section for the lower (left) and higher (right) mass models. 
Note the different scales for the contour plots. The points on the grid are shown with open squares, with the contours determined via 
interpolation. 



connection between all three quantities for both mass scales 
(unbiased cross-section only). 

First, the relation between the Einstein radius and the 
dark matter fraction shown in the bottom panels is per- 
fectly one-to-one. This fact follows from our variation of the 
dark matter profile while fixing the stellar component (see 
Appendix[B]for a derivation of this result). Second, the rela- 
tion between the lensing cross-section and the dark matter 
fraction, shown in the top panels, is not one-to-one, though 
there clearly is a mean relation with some scatter around it. 
The Einstein radius is computed in the lens plane, while the 
cross-section is computed in the source plane. The mapping 
between the two planes depends on the density profile in a 
different way than the dark matter fraction depends on the 
density profile. Therefore, variations in Cd. m and 7d m affect 
the mapping, and thus a to t, in a different way than they 
affect /d m , resulting in the observed scatter. 

While the mass range under consideration is somewhat 
different, we can compare our higher mass m odel qualita- 
tively against the results in lFedeli et al.l l|2007l ) for the effect 
of halo concentration on the strong lensing cross-section. In 
that paper, they use single-component NFW halos from sim- 
ulations to model lensing by galaxy clusters (for which this 
single-component model is not as bad an approximation as 
it is for galaxy-scale halos, where the baryonic component 
must be incorporated to obtain even remotely plausible lens- 
ing prop erties). Our results are qualitatively similar to the 
results in lFedeli et al.l (|2007T ). who find that: (a) for a sample 
with a broad mass range, the lensing cross-section is domi- 
nated by those of the highest mass and therefore lowest con- 
centration; and (b) for a sample with a narrow mass range, 
it is dominated by those clusters at the higher end of the 
lognormal scatter in halo concentration at fixed mass. Their 
results are consistent with our finding that the mass increase 
going from lower to higher mass scale has a far stronger effect 
on atot than the decrease in concentration, but that at fixed 
mass, changes in concentration can alter the lensing prob- 
abilities by ~ 50 per cent or more. Thus, the abundance 
of high-concentration halos as lenses should not be taken 
as a reflection on the underlying distribution of concentra- 
tions; and likewise, claims that too many high-concentration 



clusters have been observed should be revisited with a con- 
centration distribution from simulations modified by <T to t as 
a function of mass and concentration. 

W e also compare against the results of iKuhlen et al.l 
(2004), who use spherical, single-component NFW model 
to understand the lensing properties of cluster DM halos. 
That paper finds results similar to ours: first, that increases 
in concentration increase a to t', second, that these increases 
correlate with the mean of the image separation distribu- 
tion (and therefore -Rein), so the distribution of concentra- 
tion parameters will impact the observed image separation 
distribution for lensing systems. We can compare our Fig. 1151 
(lower left) against their fig. 6, which shows the median of 
the concentration distribution shifting by approximately 70 
per cent when accounting for lensing selection bias. Our shift 
is significantly lower than that (14 per cent) because their 
analysis uses NFW halos only whereas ours includes a non- 
negligible stellar component that will tend to reduce the 
effects of changing the NFW concentration. The different 
mass range may also play a role in the difference. 

These results have impli cations for the results in 
IComerford fc Nataraianl l|2007l ). which includes a com- 
pilation of inferred NFW profile concentrations for 
strong lensing-selected clusters and a fit to a power-law 
concentration-mass relation. Ignoring issues of possible mod- 
eling bias (for example, due to overly simple modeling of 
cluster shape or substructure), we can attempt to explain 
the fact that the normalisation of their concentration-mass 
relation is at least 20 per cent higher than in iV-body sim- 
ulations. While they consider the possibility that this high 
normalisation may be due to adiabatic contraction, we note 
that for our group-scale model (as in Fig. I15[) the median 
of the concentration distribution at fixed halo mass was in- 
creased due to selection bias by nearly that amount. There 
is of course the caveat that their mass range is some what 
higher than ours, but the results of lFedeli et al.l (|2007l ) also 
support this possibility. Consequently, adiabatic contraction 
(which has not been proven to occur in cluster-mass halos) 
may not be needed, and selection bias may be significant 
enough to explain nearly all of this apparent discrepancy 
with simulations. 
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Lower (left) and higher (right) mass spherical models 




Figure 8. For lower (left column) and higher (right column) mass scales, we show: (bottom) — The Einstein radius as a function of the 
projected dark matter fraction (within the Einstein radius) for all 80 models on the 7dm _c dm grid; and (top) — The cross-section as a 
function of the dark matter fraction. The different point colours and types are the values at fixed 7 dm : hexagons in {black, red, blue, green, 
cyan, yellow, magenta} for 7 dm = {0.00,0.54,0.88,1.14,1.35,1.52,1.66}, and crosses in {black, red, blue} for 7 dm = {1.79,1.90,2.00}. 
The symbol sizes from smallest to largest indicate the 8 increasing values of e dm at fixed 7d m - 



5.2 Varying concentration in the triaxial case 

We explore the variation of orientation-averaged cross- 
sections with dark matter concentration more coarsely for 
triaxial models, to verify that the trends remain roughly the 
same as for spherical models. This assumption is acceptable 
for our coarse grid in 7d m as presented in Section 14.11 but 
we need to confirm this for Cdm as well. Because of the need 
for 24 viewing angles for any triaxial model, we use only one 
value of 7d m = 1, and we restrict the grid in Cd m to five 
values ranging from —2a to +2a in steps of la. Thus for 
the lower mass scale, we go from Cdm = 4.2 to Cdm = 16.8; 
for the higher mass scale, from Cdm = 2.8 to Cd m = 11.2. 
The actual values at the fiducial concentration are given in 



Table [2] for all magnification bias modes, so we only present 
trends relative to those values. 

There are several points we infer from the resulting 
Fig. [9] First, our conclusions from the spherical case are, 
broadly speaking, robust for the triaxial case: variation of 
Cdm with all other parameters fixed leads to enhancements 
of the total lensing cross-section by factors of a few. In par- 
ticular, the enhancement can be factors of ~ 2 for the lower 
mass scale, and ~ 4 for the higher mass scale, where the 
dark matter component is more important. The reason the 
range of enhancements is smaller than in Section [S7T] is that 
the range in Cdm is smaller; however, the range of enhance- 
ments is in line with the spherical case when we restrict the 
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Figure 9. For lower (left column) and higher (right column) tri- 
axial mass models, we show the variation of the total, orientation- 
averaged lensing cross-section (normalised to the value at the fidu- 
cial concentration) with concentration in the top panels, and the 
variation of the 4:2 image ratio with concentration in the bottom 
panels. Line colours and types indicate the magnification bias 
mode, as described in the caption for Fig. [l] 



in which lensing can occur, but also to the scaling of the 
surface mass density with radius within this region. 

Despite some subtle differences between changing con- 
centration and inner slope, it is clear that in strong lensing- 
selected analyses, we expect significant selection biases (fac- 
tors of several) favouring steeper inner slopes 7d m and higher 
concentrations Cd m at fixed mass. This selection bias is man- 
ifested not just in a to t but also in 0-4/0-2. This result im- 
plies that an analysis that recovers the distribution of inner 
slope and/or concentration with high statistical precision 
does not measure the intrinsic distribution of these quanti- 
ties directly, but must account for the selection biases first 
using simulations. 



6 DEPROJECTED SERSIC DENSITY 

In this section, we investigate selection bias with density pro- 
file parameters for the composite deprojected Sersic density 
profiles. Previous work has characterised the basic lensing 
prop ert i es of spherical, s i ngle-component d eproj ected Sersic 
models l|Cardonel 12004 iBaltz et all l200Sh . but we study 
the full lens statistics of two-component spherical and non- 
spherical models. We divide this analysis into three parts, 
where we first address the comparison against cusped mod- 
els in the spherical case; we next consider orientation depen- 
dence in the mixed and triaxial cases; and we finally allow 
the Sersic index to vary separately for the dark matter and 
stellar components. 



grid for the calculations with spherical models in the same 
way. 

One point we were not able to test using the spheri- 
cal case is the issue of the quad/double ratio 0-4/0-2- Here, 
we see that this number also can change by factors of a few 
for reasonable ranges of concentrations for both mass scales. 
This change presumably results from the fact that the dark 
matter is intrinsically less round than the stellar compo- 
nent, so increasing the dark matter concentration increases 
its importance in the inner parts and therefore increases the 
projected ellipticity of the total mass distribution. 

We can also contrast Fig. [9] with the "triaxial" panels 
of Figs [3] and 3] showing the variation of respectively o to t 
and 04/02 with 7dm at the fiducial concentration. First, the 
increase in total lensing cross-section when going from the 
median concentration to that +2a is about half the increase 
when going from 7dm = 1 to 7d m = 1.5. Second, there are 
some subtleties that mean that increasing 7d m and increas- 
ing Cdm do not affect the lensing cross-section in quite the 
same way. In particular, the ranking of magnification bias 
modes is different in the two cases: when increasing 7dm, the 
otot for the higher flux limit magnification bias modes does 
not increase as much as for the unbiased case, whereas the 
opposite is true when increasing Cdm- Furthermore, 0-4/02 
tends to flatten out when going from 7 dm = 1 to 7 dm = 1-5, 
whereas it increases by a factor of ~ 2 when going from the 
median concentration to +2o. Both of these differences be- 
tween increasing 7d m and increasing concentration are due 
to the fact that the biased cross-sections and the four-image 
cross-section are related not just to the size of the region 



6.1 Comparison against cusped models 

We first ask whether deprojected Sersic models lead to 
significantly different lensing cross-sections compared to 
cusped models, for both mass scales in the spherical case. 
This question is important for example in theoretical calcu- 
lations of the expected number of lenses from a particular 
strong lensing survey. 

The fiducial cusped models with 7d m = {0.5, 1.0, 1.5} 
have projected dark matter fractions within the correspond- 
ing Einstein radii of /d m = {0.11,0.27,0.62} (lower mass) 
and /dm = {0.18, 0.49, 0.88} (higher mass). The fiducial de- 
projected Sersic model has / dm = 0.55 and / dm = 0.85 
for lower and higher mass scales, respectively, which means 
that qualitatively we expect the Einstein radii and cross- 
sections to be more comparable to the 7 dm = 1.5 cusped 
model than to the 7dm = 1 case. The results from analytic 
GRAVLENS calculations show that this expectation is correct: 
for the lower mass scale, otot = 1.57arcsec 2 , comparable to 
otot = 1.40arcsec 2 for the 7dm = 1.5 case, while o to t = 
0.53 arcsec 2 for 7d m = 1. Likewise, for the higher mass scale, 
otot = 10.5 arcsec 2 , which is closer to o to t = 15.9 arcsec 2 
for the 7dm = 1.5 case than to o to t = 1.63 arcsec 2 for the 
7dm = 1 case. All numbers quoted here are unbiased cross- 
sections, and the trends for i? c i n are similar. We conclude 
that for the same mass within spheres of over-density 180/9, 
deprojected Sersic models with observationally-motivated 
Sersic indices and half-mass radii tend to give higher lens- 
ing cross-sections than cusped models with the "canonical" 
inner slope value 7d m = 1. It is clear that the selection bias 
with mass for the deprojected Sersic models is broadly con- 
sistent with what one would expect for the cusped models 
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with 7d m = 1.5, which is stronger than the selection bias for 
the 7dm = 1 models. 

When we include the effects of magnification bias, we 
find that the enhancement in the cross-section is similar to 
the enhancement for the 7dm = 1-5 cusped models. The 
image separation distributions are fairly narrow as for the 
cusped models, with a mean of 1.98 iicin (lower mass) and 
1.88 i? c in (higher mass). This result is consistent with our 
finding in Paper I that the composite deprojected Sersic 
model for the galaxy scale is very close to isothermal, 
whereas the higher mass scale is more dominated by the 
dark matter and therefore is somewhat less isothermal. 



6.2 Orientation-dependence 

We had found, in Section f3. II that for the triaxial and mixed 
shape cusped models, the total unbiased cross-section <T to t 
can vary by factors of several when we change the viewing 
angles. Likewise, the ratio 0-4/0-2 can vary from near zero to 
~ 15 per cent. 

We now compare this result for the orientation depen- 
dence in cusped models against the results for the depro- 
jected Sersic models. We find that the results are very sim- 
ilar for variation of both o to t and of 04/02, which suggests 
that these numbers do not depend strongly on the func- 
tional form of the density profile (for reasonable models). 
Furthermore, when we include the effects of magnification 
bias, we find that it changes o to t and 0-4/0-2 in similar ways 
as it does for the cusped, 7<jm =1.5 models. This change 
includes sharp increases in o tot and significant enhancement 
of the quad/double ratio. 

When averaging over orientation, we also find that, as 
for cusped models, the value of <7 tot does not depend strongly 
on halo shape (when comparing triaxial, mixed, and spher- 
ical models). 

6.3 Sersic index 

In this section, we consider varying the Sersic indices of both 
the dark matter (nd m ) and stellar component (n*), in the 
spherical case, for both mass scales. All variation is done 
at fixed mass. Because an increase in Sersic index leads to 
a more highly concentrated profile in the inner regions, we 
expect that with appropriate normalisation (discussed in Pa- 
per I), the lensing cross-section will be an increasing function 
of Sersic index. The questions we address are (1) the extent 
of the selection bias that is possible for a reasonable range of 
Sersic indices, and (2) the relative importance of the Sersic 
index for each component of the profile. We adopt the ranges 
in ridm and n* as derived in Paper I, and sample five values 
equally in both log and log n& m , resulting in a 5 x 5 grid. 
For ridm and n* values different from the fiducial values (in 
the centre of the grid), we change po to maintain the total 
mass and scale radius. 

Fig.[lO]shows the unbiased cross-section as a function of 
the Sersic indices for both mass scales, with contours in the 
base- 10 logarithm relative to the cross-section at the fidu- 
cial grid point, i.e., log (crtot/otot.fid). The increase in cross- 
section with both ndm and n* is physically understandable. 
As shown, the total variation on the grid is from ~ 0.5 to 
2ctot,fid for the lower mass scale (factor of 4), or ~ 0.45 
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Figure 10. Contour plots of the base-10 logarithm of the unbi- 
ased cross-section for the lower (top) and higher (bottom) mass 
scales, relative to the cross-section at the fiducial grid point. 

to 2.6 otot, fid for the higher mass scale (factor of 6). While 
these increases are significant, they are more comparable 
to the changes in cross-section for the cusped models when 
varying Cd m than varying 7d m - The reason for this is that 
changing 7d m for a cusped model causes a large change in the 
inner parts of the intrinsic and projected profile, but chang- 
ing the Sersic index leads to substantial changes in both the 
inner and the outer intrinsic profile that partially cancel out 
when projected along the line-of-sight. The shape of the con- 
tours, which run diagonally at low rid m and n*, but closer 
to horizontally at high rid m and vertically at high n*, is also 
physically understandable: at low ridm and n*, both compo- 
nents are roughly equally important, whereas at high nd m , 
the stellar component (and therefore n*) is less important, 
and at high n*, the dark matter component (and therefore 
ndm) is less important. The effect at high 7id m is more ex- 
aggerated for the higher mass scale, because this model was 
designed to be more dark-matter dominated than the lower 
mass scale model. 

Fig. llll shows the trends of Rdn/R'h.-t f° r the lower mass 
scale on this grid. Note that since we normalise the profiles 
by fixing the scale radius and mass while varying po, the 
scale radius R' h i < is the same at each point on the grid, so 
only R e i n is varying. The exceptionally horizontal contours 
require some explanation, since the cross-section contours 
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Figure 11. Contour plot of Rdn/Rh * for the lower mass dcpro- 
jected Sersic models. 

in Fig. [TU] are not strictly horizontal. This result can be ex- 
plained by the deflection curves, shown in Fig. [12] for this 
mass scale. In this plot, the Einstein radius is where the 
one-to-one line crosses the a(R') curves for the composite 
models, and as we have already seen, these crossings occur 
in nearly the same place. This results from our normalisa- 
tion convention, which fixes the same half-mass radius R' h + 
for each of the n* values, meaning that it fixes the average 
convergence K(R' h , i ,) and therefore the deflection at R'h,*, 
a = R' h it K(R' h to the same value. Because the DM com- 
ponent is also fixed, it means that the models with different 
n* have the same R B1U and the same f' Aul within i? c i n . How- 
ever, this is purely a consequence of our chosen redshifts and 
the balance of stellar and dark matter in our models; had we 
arranged it so that i?em <C R'h,*, then n* would indeed af- 
fect i?cm- However, the cross-sections do not have to be the 
same despite identical R cln values for varying n*, because 
they are set by the outer (source-plane) caustic, which is re- 
lated to the inner (lens-plane) critical curve, determined by 
the slope of the deflection on small scales. Clearly, in Fig. 1121 
while the deflections are the same at R c in, their slopes are 
quite different within 7? e in, leading to different lensing cross- 
sections. Finally, we note that despite our superimposing two 
non-trivially complicated deprojected Sersic models to make 
the composite galaxy models, the deflection curve for the 
composite model is remarkably close to the flat, isothermal 
case (a = .Rem), which we have also seen with the cusped 
models. This fact explains the very narrow image separation 
distribution for all of these composite models. 

Finally, we consider the selection bias trends for biased 
cross-sections, particularly the one using the highest flux 
limit and total magnification. That trend is presented for the 
lower mass scale in Fig. 1131 and shows that the shape of the 
contours is quite different from the unbiased case (Fig. 1101 
top panel). While Udm has the same, significant effect on 
the cross-section as in the unbiased case, the effect is much 
less significant for n* and has changed sign: stellar compo- 
nents that are less concentrated give slightly higher biased 
cross-sections. While there is no simple, physical explana- 
tion for this result, it is apparently related to the scaling of 
the convergence within the Einstein radius, and the fact that 
a very concentrated (high n*) stellar component only has a 
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Figure 12. Deflection as a function of scale for the lower mass 
deprojected Sersic mass model. The black dashed line is for the 
fiducial DM component; the black, red, and blue dotted lines are 
for the stellar components with the lowest, fiducial, and highest 
n*; and the black, red, and blue solid lines are the correspond- 
ing composite models (stellar plus DM). Finally, the green solid 
line corresponds to a = R'; the intersection of this line with the 
deflection curves occurs at R eln . 
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Figure 13. Contour plot of the base-10 log of the biased cross- 
section for the lower mass scale, relative to this cross-section at 
the fiducial grid point. 



very small region where the convergence is extremely high, 
whereas a less concentrated one has a larger area where it is 
fairly high, which turns out to be more important for these 
very biased cross-sections. 
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7 ROBUSTNESS OF CONCLUSIONS TO 
CHANGES IN REDSHIFT 

Finally, to ascertain how much our conclusions about selec- 
tion biases in the previous sections depend on our chosen 
fiducial lens and source redshifts, we now try varying the 
redshifts in several ways. 

If we vary both the lens and the source redshifts in a 
way that preserves E c , then the physical scale within the 
galaxy that is measured by strong lensing is the same, but 
it will lead to different image separation distributions, scaled 
by the angular diameter distance Dl of the lens galaxy. This 
change may lead to observational selection effects that de- 
pend sensitively on the survey image search strategy (res- 
olution and maximum angular separation), but no physical 
selection effects, so we do not attempt to model this in any 
general way. 

Now consider fixing the lens redshift (and therefore the 
conversion from physical to angular scale), but shifting the 
source redshift to rescale the critical surface mass density 
E c . This shift will change the Einstein radius R c i n and there- 
fore the physical scale of the density profile that is probed 
by strong lensing. Thus, this change may result in different 
physical selection biases. 

If we fix the lens redshift to our fiducial zl = 0.3 but 
move the source redshift considerably lower, from the fidu- 
cial zs = 2.0 to z' s = 0.6, then the critical surface density in- 
creases by ~ 70 per cent. This change means that the strong 
lensing is due to smaller scales, where the average projected 
surface density is higher, and that the lensing cross-section 
for this density profile will be lower. Because of the impor- 
tance of smaller scales, the projected dark matter fraction 
/dm wm decrease. As a test, we recompute R c i n and the 
lensing cross-sections for both mass scales as a function of 
the dark matter halo inner slope 7d m to see the degree to 
which the selection bias we identified with our fiducial lens 
and source redshifts persist. We also increase the resolution 
of the surface density maps used by GRAVLENS for these cal- 
culations by a factor of two, while maintaining the box size, 
to ensure that the critical curves can be properly resolved. 

We also consider a value of E c that is 40 per cent lower 
than our fiducial value. While this low E c does not corre- 
spond to any sensible value of zs for our fiducial zl, we 
assume a higher zl ~ 0.6 and zs ~ 3. Because of the dif- 
ferent zl, all angular scales output by our pipeline must be 
rescaled by D A {z = 0.3)/ D A (z = 0.6) = 0.66, i.e., the Ein- 
stein radii must be rescaled by 0.66 and the cross-sections 
by 0.66 2 . However, the relevant quantities for selection bias 
associated with true variation of galaxy density profiles, e.g. 
/dmi are preserved. Since we expect increased Einstein radii 
and cross-sections, we preserve the resolution but increase 
the box size by a factor of 2. 

Table [3] shows what happens to relevant quantities for 
lensing when we vary E c to values 70 per cent higher and 40 
per cent lower. We consider only the unbiased cross-sections 
for the spherical shape model. First, we can see that when 
we increase (decrease) E c , we decrease (increase) R cln and 
/d m , as expected. The change in 7? e m can be significant, 
up to a factor of ~ 2 for the changes in E c that we have 
tested. We find that the selection bias with 7 dm , represented 
in the last two rows of this table, is nearly independent of 
E c for 7dm below the fiducial value of I, and for higher 7d m 
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Figure 14. Critical surface density, plotted as 
log [£ c /(M0/pc 2 )], as a function of lens and source red- 
shift. For zs < zl, where S c is formally infinite, we have set it 
to the maximum value shown, 10 4 Mq/pc 2 . Asterisks indicate 
the minimum source redshift for a given lens redshift for which 
we have confirmed that our conclusions about selection bias are 
valid. 



it is a slightly increasing function of E c . The sign of this 
change is surprising, in the sense that if f' Am is lower than 
we might naively expect the parameters of the dark matter 
distribution to be less significant. However, we can explain 
this change by the fact that the cross-section a tot is deter- 
mined by the outer (source-plane) caustic, which relates to 
the inner (lens-plane) critical curve. This in turn is deter- 
mined by the slope of the density profile on small scales, so 
if the lensing properties are determined by smaller (larger) 
physical scales, as when we increase (decrease) E C) then 7d m 
becomes more (less) important and selection bias becomes 
stronger (weaker). 

However, it is encouraging that the selection bias only 
changes by +14 and —4 per cent (+41 and —13 per cent) 
for changes in E c of +70 and —40 per cent for the lower 
(higher) mass model. This result suggests that within the 
range of E c that we have considered in this subsection, the 
conclusions of our paper regarding physical selection bias 
are not strongly dependent on our choice of fiducial lens 
and source redshifts. (The same is not true for observational 
selection biases, which may be strongly affected by changing 
zl and thus the characteristic angular scale of the system.) 
On Fig. 1141 we have indicated the lens and source redshift 
combinations for which our conclusions are robust according 
to the tests in this subsection. 

We must also consider the issue of magnification bias. 
Magnification bias complicates a change in the redshift of 
the system in two ways. First, if the source redshift is dif- 
ferent, then for a given survey flux limit, the luminosity 
threshold is different. For a broken power-law luminosity 
function that is steeper at the bright end, this means that if 
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Table 3. Summary of the lensing results for the cusped, spherical model when varying S c so that different physical scales are probed. 
The factor of 0.66 in R e in for the lower E c cases is shown separately since it is not indicative of the physical scale within the galaxy. 



Quantity Fiducial S c Higher S c Lower E c 





Lower 


mass scale 








-Rein (7dm = 0.5) 




0.58 


0.36 


0.79 x 0.66 


= 0.52 


Rein (7dm = 1-0) (arcsec) 




0.68 


0.41 


0.95 x 0.66 


= 0.63 


-Rein(7dm = 1.5) 




1.04 


0.61 


1.46 x 0.66 


= 0.97 


f'dm( R ' <Rein,7dm=0.5) 




0.13 


0.08 




0.17 


fdm( R ' < R ein, 7dm = 1.0) 




0.28 


0.20 




0.37 


/ dm (R' < Rein, 7dm = 1.5) 




0.60 


0.51 




0.67 


O"tot(7dm = 0.5)/(Ttot(7dm = 


1.0) 


0.79 


0.79 




0.79 


ftot(7dm = l-5)/<xtot(7dm = 


1.0) 


2.65 


2.96 




2.51 




Higher 


mass scale 








Rein (7dm = 0.5) 




1.18 


0.63 


1.73 x 0.66 


= 1.14 


Rein (7dm = 1.0) (arcsec) 




1.75 


0.88 


2.75 x 0.66 


= 1.82 


Rein (7dm =1.5) 




4.44 


2.33 


6.75 x 0.66 


= 4.46 


fdm( R ' < Rein, 7dm =0.5) 




0.22 


0.13 




0.30 


C(K' <Rein,7dm = l-0) 




0.51 


0.36 




0.63 


/dm( R ' <Rein,7dm = l-5) 




0.87 


0.80 




0.91 


0"tot(7dm = 0.5)/CTtot(7dm = 


1.0) 


0.52 


0.51 




0.52 


0"tot(7dm = l-5)/(Ttot(7dm = 


1.0) 


9.19 


12.9 




8.04 



the source is more distant, the limiting luminosity becomes 
brighter, so the actual number of lensing systems is mod- 
ulated due to both the different number density of sources 
above that limit and by the change in effective a to t due to 
the steeper luminosity function at the limit. Second, since 
the source luminosity function may be a function of red- 
shift, we anticipate that the biased cross-sections, which are 
the relevant quantity for any particular survey given its flux 
limit and search strategy, may also change due to this effect. 
For our analysis in this section, we have ignored this change 
in the magnification bias formalism, but to determine the 
cross-sections for a real survey this evolution would have 
to be taken into account. In any case, both of these mag- 
nification bias effects qualify as observational rather than 
physical selection biases that are more properly simulated 
in the context of a particular strong lensing survey. 



8 CONCLUSIONS 

Using the pipeline developed in Paper I, we have presented 
a series of tests for selection biases in point-source strong 
lensing surveys due to physical differences between galaxy 
or group mass systems. These tests reveal whether the distri- 
bution of physical parameters, such as profiles and shapes of 
(projected) densities, in strong lensing systems reflects that 
distribution among the galaxy population in general. 

We have attempted to distinguish between physical se- 
lection biases which relate to the parameters of the galaxy 
density, versus observational selection biases which depend 
on the lensing detection strategy (specifically, the flux limit, 
angular resolution, and maximum angular separation for de- 
tection of a lensing system). The former can be understood 
for lensing surveys in general, whereas the latter need to be 
simulated on a case-by-case basis. 

To organize our conclusions, we first give a sample ap- 
plication to illustrate how our analysis may be applied to 
real lens surveys, then give a more complete list of our con- 



clusions about selection biases in strong lensing surveys, and 
finally give a more general discussion of our results. 



8.1 Example application 

As an example of how this work may be applied to the anal- 
ysis of real data, we consider the joint distribution of Cd m 
and 7dm and show how the distribution that would be ob- 
served in some strong lensing survey is related to the in- 
trinsic distribution. In this example we use the spherical 
higher mass model. The concentration is known from N- 
body simu lations to have a lo gnormal distribution at fixed 
mass (e.g.. lBullock et ai]|200ll ). The distribution of the loga- 
rithmic slope is not well known, however, because of difficul- 
ties in resolving small scales. Indeed, some recent work has 
suggested that DM halo profiles do n ot converge to a sin - 
gle asymptotic logarithmic slope (e.g.. iNavarro et al.ll2004f ). 
although even in that case one can imagine constraining 
(in both simulations and real data) the inner slope at some 
fixed fraction of the virial radius and calling that 7dm- We 
assume, for illustration, that this procedure may give a mean 
(7dm) = 1 and a Gaussian scatter of 0.25. We expect that 
any constraints on 7d m and Cdm based on measuring en- 
closed mass will cause 7dm and Cdm to be anti-correlated, so 
we assume a correlation coefficient of —0.3 between 7d m and 
ln(cdm). We emphasize that the 7d m distribution and the 
correlation coefficient are merely guesses made for the pur- 
pose of this exercise, which means that our results should not 
be taken as a prediction for what would be seen in reality; 
they are intended just to be a useful illustration. 

Fig. [15] shows our assumed intrinsic distribution (lower 
right panel), along with the distribution that would be ob- 
served among lens galaxies (upper right panel). The two 
distributions are related via the lensing selection function 
shown in Fig. [7] (For this exercise we use the selection func- 
tion without magnification bias.) Lensing selection effects 
clearly shift the peak of the joint distribution to higher val- 
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Figure 15. Illustration of the procedure for removing the effects of lensing selection bias from the observed distribution of parameters, 
shown for the spherical, higher mass model. The top left panel shows the joint distribution psL(7dm> Cdm) a s it may be observed in some 
strong lensing survey. The bottom right panel shows the corresponding intrinsic distribution Pint (7dm > c dm)i which would be recovered 
from the observed distribution by removing the selection bias function shown in Fig. [7] In both cases the grayscale and contours are 
logarithmic, and the horizontal and vertical lines are drawn at the median values of 7 dm and c dm for the intrinsic distribution. We 
emphasize that the instrinsic joint distribution is not well known from observations or simulations, so we have constructed a simple 
example for illustration (see text for more discussion); thus, this figure should not be taken as a prediction for what will be seen in 
lensing surveys. The top right and bottom left panels show cumulative distributions for 7 dm and Cd m (respectively) obtained by taking 
slices through the joint distributions the intrinsic median values of Cdm or 7dm- In each case the dashed line represents the "observed" 
distribution while the solid line represents the recovered intrinsic distribution. The dotted horizontal lines indicate the median and ±lcr 
range. 



ues of 7dm and c dm . To further quantify the effects, we also 
show 1-d slices through the joint distribution at fixed val- 
ues of 7dm or c dm - These plots demonstrate that the 7d m 
and Cdm distributions are both shifted and broadened in 
the "observed" distribution compared with the intrinsic dis- 
tribution. In particular, the mean and median values are 
shifted by 24 (7dm) and 14 (cdm) per cent. Consequently, 
any attempt to recover the mean of Cdm and/or 7dm for 
some sample of strong lenses need to be corrected for se- 
lection bias before comparing against predictions that were 
made for a more general galaxy sample. 



8.2 Summary 

The biases in 7d m and Cd m are just one example of the phys- 
ical selection effects we have explored. Here is a more com- 
plete list of our main conclusions: 

• There is a strong preference for lensing systems to be at 
higher mass. This trend is counterbalanced by the strongly 
decreasing mass function at high mass, but will nonethe- 
less significantly increase the mean mass for lensing systems 
relative to that for all galaxies. The degree of preference de- 
pends sensitively on how the mass is distributed in the inner 
parts of the galaxy, and on the details of the magnification 
bias (flux limit and detection strategy). [Section 12. 4j 

• Lensing systems show a significant preference for cer- 
tain orientations, with factors of 2-3 variation in lensing 
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cross-sections depending on their shape. In particular, the 
total cross-section is maximised when the surface mass den- 
sity is maximised, i.e., when the galaxy is viewed with the 
long axis along the line of sight; and the four-image cross- 
section is maximised when the projected ellipticity is max- 
imised, i.e., when the intermediate axis aligns with the line 
of sight. [Sections 13. II and 16. 2] 

• The total, orientation-averaged lensing cross-section is 
not strongly dependent on the intrinsic shape of the galaxy. 
In contrast, the quad/double ratio may be a powerful probe 
of the distribution of the shape of the mass density in the 
inner parts of galaxies, although it also has a strong depen- 
dence on observational issues such as the source luminosity 
function. [Section |3. 2] 

• An observational scarcity of three-image lens systems 
may place a strong constraint on the abundance of strongly 
oblate matter distributions. [Section 13. 2| 

• Statistical variations in the dark matter density proper- 
ties, such as the inner slope 7<jm and the concentration Cdm, 
may lead to a strong selection bias amounting to factors of 
many in a to t, and to tens of per cent variation in 04/(72, 
nearly independent of the shape of the galaxy model and 
the observational issues relating to the magnification bias. 
Consequently, even if analyses are able to break the degen- 
eracy between these parameters, for example using strong 
lensing plus some other technique to measure the profile at 
larger scales, they still cannot claim to measure the intrinsic 
distribution of these quantities directly. [Sections 14.11 15.11 
and EH 

• Adiabatic contraction has a less significant effect on the 
lensing cross-section than increasing the dark matter inner 
slope for a generalised NFW model from 7d m = 1 to 7d m 

I Si However, the increase may still be a factor of two, which 
means that if galaxy formation scenarios dictate whether AC 
occurs, there may be a selection bias towards those galaxies 
that have undergone significant adiabatic contraction, and 
against those that have not. [Sections 14. II and 14. 3] 

• The lensing cross-sections of deprojected Sersic mod- 
els depend on the indices n* of the stellar and rtd m of the 
dark matter component. Higher Sersic indices correspond to 
a more concentrated profile and therefore a higher lensing 
cross-section at fixed mass. Reasonable changes of these in- 
dices, within the observational scatter, can change the lens- 
ing cross section by factors of several. For the fiducial Sersic 
model, 7? e in and atot are higher than for cusped models with 
the "canonical" inner slope value 7d m = 1, and closer to 
those for the 7dm = 1.5 cusped model. [Section [6] 

• While varying the lens and source redshifts changes the 
range of physical scales probed using strong lensing, we have 
found that the dependence on E c is not very strong, and 
therefore that our conclusions regarding the extent of selec- 
tion biases are valid for a wide range of lens-source redshift 
combinations (Fig. I14|l . Furthermore, our results suggest 
that the physical selection bias is modified (when changing 
lens and/or source redshift) such that the very inner part of 
the profile is more (less) important if S c is higher (lower), 
making selection bias more (less) important than for our 
fiducial combination of lens and source redshifts. [Section [7] 



8.3 Discussion 

Where possible, we have compared our conclusions on phys- 
ical selection biases with earlier papers in the relevant sec- 
tions. In general our conclusions are consistent with those 
of previous studies, but our pipeline allows a more uni- 
fied approach to selection biases in lensing surveys while 
using realistic galaxy models. This unified approach al- 
lows for non-trivial extensions of previous results, to two- 
component models, different profile types, and more shapes 
than were used in previous studies. Furthermore, we can 
understand trends seen in strong lens analyses based on 
N-body plus sem i -analytic modeling (|Hilbert et al.l 120071 ; 
IMoller et aT]|2007l ; iHilbert et afll2008l ) in terms of the inter- 
play between different t ypes of selection bia ses. For example, 
the tendency found by IMoller et al.l (|2007h for quad lenses 
to have a pronounced disk component comes from an inter- 
play between mass bias, shape bias, and orientation bias, 
as discussed in Section 13.2.21 A comparison with previous 
results also teaches us lessons about the limitations of these 
simple models: for example, that when estimating selection 
biases with dark matter halo parameters, we must include 
a baryonic component so that we do not overestimate the 
severity of the selection bias, as discussed in Sections 14.11 
andO 

Finally, our work has revealed a number of useful points 
related to the interpretation of observational strong lensing 
results in general. For example, despite the relative complex- 
ity of our mass models, the image separation distribution in 
all cases has a mean near 2i? c i n and a fairly small width 
(Section 14. 2|) . as in the isothermal case. The dependence of 
this distribution on magnification bias is a very small effect. 
On a related note, these models are very close to isother- 
mal near i? c i n despite strong variations at smaller or larger 
scales (see paper I figure 5), so a measurement to constrain 
the slope of the projected density profile at R e i n does not 
allow discrimination between profiles with different 7d m or 
Cdm, unless another technique is used to constrain the pro- 
file independently on much smaller or larger scales (see also 
Paper I). Note, however, that because of the significant devi- 
ations from isothermality on other scales, these models have 
very different relationships between R c [ n and fftot than for 
an isothermal model (see also Table [TJ. 

We can apply several of o ur conclusions to previously- 
published work. For example, iKoopmans et all (|2006l ) use 
a joint lensing and kinematics analysis of 15 SLACS lenses 
to infer that the average slope of the 3d density profile near 
the Einstein radius is (7) = 2.011q'q3 (68 per cent CL), with 
intrinsic rms scatter of 0.12. We defer the consideration of 
possible modeling biases to future work, and consider the 
possible impact of selection biases here. This extremely high 
level of isothermality is, at first glance, surprising; but then 
we realise from paper I figures 5 and 7 that in fact many 
different combinations of stellar and DM models can con- 
spire to yield a total model that is nearly isothermal over a 
range of scales including the Einstein radius. Consequently, 
a measurement of the total density profile in the vicinity of 
the Einstein radius cannot, on its own, u niquely determine 
the fo rm of the DM profile. However, once IKoopmans et al.l 
(2006) impose a prior on the stellar mass-to-light ratio they 
are able to constrain the projected dark matter fraction, 
which (in combination with knowledge of the projected lu- 
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minosity profile) might provide constraints on the inner dark 
matter profile. However, due to the partial degeneracy be- 
tween mass, 7d m and Cd m in strong lensing properties, this 
constraint does not allow for an explicit determination of 
the dark matter profile parameters without additional con- 
straints from larger radii. In other words, there are many 
models in that three-dimensional space that would give the 
observed J? m and / dm . 

There are a number of selection bias-related problems 
that involve a higher level of complexity and that we have 
therefore reserved for future work. These include, but are not 
limited to: (i) the effect of varying stellar mass-to-light ratios 
for a given galaxy model; (ii) modeling satellite galaxies that 
have their own stellar and dark matter components but are 
distributed within the halo of a group or cluster (rather than 
residing at the centre), and may have the dark matter and/or 
stars tidally stripped; (iii) lensing of extended sources rather 
than point sources, for which the basic trends should be 
qualitatively similar, but for which observational selection 
effects and observing conditions are crucial; (iv) more use- 
ful ways of quantifying image separation distributions for 
systems with more than two images, and what information 
can be learned in this way; (v) intrinsic misalignment of the 
non-spherical dark matter and stellar components, warps, 
and other more complex shape models. 

In addition to this future work on selection biases, we 
also have designed this pipeline so that it can be easily 
adapted to study modeling biases, which can occur when us- 
ing overly simplified mass models. Work on modeling biases 
can be extended to include strong + weak lensing analyses, 
as can be done at cluster mass scales. With the addition of 
kinematics, we can study modeling biases not just in strong 
lensing systems, but also in kinematic and in kinematic + 
strong lensing analyses. In short, this pipeline should be a 
useful tool in understanding the analysis requirements for 
current and future strong lensing, kinematics, and weak lens- 
ing surveys that aim to understand galaxy density profiles. 

Finally, we note that considerable additional work 
would be necessary to use these results to estimate selec- 
tion biases in a real survey. For example, observational se- 
lection biases need to be accounted for. This may include 
the imposition of a maximum or a minimum angular separa- 
tion between images depending on the use of an aperture or 
a resolution-limited approach, respectively. For a n example 
of wor k considering these and related issues, see iKochanekl 
(1991). While this may change the effects of physical selec- 
tion biases, either enhancing or counteracting them, it can 
also alter the lens redshift distribution by favouring higher or 
lower redshifts, respectively, due to the conversion between 
angular and physical scale. Furthermore, the evolution of the 
survey flux limit with the source redshift distribution must 
be taken into account: we have considered a fixed luminosity 
limit, but a typical flux-limited survey would have a variable 
luminosity limit that gets brighter for higher redshift, alter- 
ing the biased cross-sections as a function of source redshift 
even if the source luminosity function is assumed to be inde- 
pendent of redshift. This will tend to push the source galaxy 
population to higher redshifts on average. 

Another point is that we have not attempted to model 
the full intrinsic distribution of galaxy properties, pint (a?)- 
We have studied instead the (physical) mapping function 
o~(x) that determines the final observed strong lensing popu- 



lation from the intrinsic one. To predict strong lensing popu- 
lations would require a full simulation of the intrinsic galaxy 
property distribution pi nt (a;), which is clearly beyond the 
scope of this paper, but which has been attempted in limited 
form using N-body s i mulations plus semi -anal ytic models by 
iHilbert et all <|2007l ). lMoller et all (|2007h . and lHilbert et all 
(2008). In short, a full simulation of the lens galaxy popu- 
lation and source quasar population would be necessary to 
estimate the parameters of the strong lensing galaxy popula- 
tion in detail for a particular survey. Nonetheless, our work 
is an important first step in unraveling and disentangling 
physical sources of selection biases for current and future 
strong lensing surveys. 
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APPENDIX A: VARIATION OF LENSING 
CROSS-SECTION WITH INNER SLOPE AND 
ORIENTATION 



In Section [4. II we quantified a selection bias based on dark 
matter halo inner slope using ratios 



{o-jv(7dm = 1-0)) 

(o-]v(7dm = 0.5)) ' 



(Al) 



where the averages indicate an average over orientations, 
and TV can indicate either the total cross-section or the cross- 
section for some number of images N. We can also define 
a comparable quantity ((i?i.s,i(iV))) for 7d m = 1.5 versus 
7dm = L Ratios of one indicate an absence of selection bias, 
and we have generally found these ratios to be > 1 (the 
strong-lensing cross-sections increase with 7dm). 

We now address the question of whether these ratios 
are the same for each viewing direction, i.e. we consider 

u c at 'i u' i \ crjv (7d m = 1.0, a'/a, 6'/ a ) /Ar>\ 

iii,o.5 (N, a /a, b a) = y / A2 

crjv(7dm = 0.5, a 1 1 a, b'/a) 

and a comparable quantity #1.5,1 (TV, a' /a, b'/a). If these ra- 
tios stay (nearly) constant while changing a'/a and b'/a, 
then the selection bias with dark matter halo inner slope is 
not a (strong) function of viewing direction for galaxies with 
non-spherical shapes. 

For the triaxial lower mass models, these ratios are in- 
deed approximately constant (to within ~ 10 per cent) for 
fftot, < J 2, and (74, independent of magnification bias mode. 
For the triaxial higher mass model, where dark matter is 
more important, we find a stronger dependence on view- 
ing direction. For 02, there are up to ~ 30 per cent level 
trends for #1,0.5 and #1.5,1 to be lower (i.e., closer to unity 
and hence indicating less selection bias with 7dm), when 
viewing the model along the intermediate axis. This change 
occurs because in that orientation, the projected ellipticity 
and therefore 0-4 is maximised, which decreases 02 because 
2-image systems can only occur outside of the tangential 
caustic. For 04, this is the case when viewing the model 
along the short axis (for #1,0.5) and along the long axis (for 
#1.5,1). 

It is not clear that the details of this orientation- 
dependence of the selection bias with 7d m is very important 
except for special situations. When analysing data from a 
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typical survey, we expect that galaxies are oriented ran- 
domly (before the imposition of strong lensing "orienta- 
tion bias"), so the angle-averaged selection biases from Sec- 
tion [47T] are the relevant quantities. We include the analysis 
here as an indication that if one isolates a special subset of 
strong lenses based on their projected shape, such as those 
for which the lens is particularly round or flattened, then 
the treatment of selection bias with 7d m and possibly other 
properties must be specialised to that subset rather than 
using the orientation-averaged case. 



APPENDIX B: RELATIONSHIP BETWEEN 
DARK MATTER FRACTION AND EINSTEIN 
RADIUS 

We seek to understand why the relation between dark matter 
fraction and Einstein radius shown in Fig. [8] is perfectly 
one-to-one. We begin with the Einstein radius definition, 
K(-Rejn) = 1j meaning that the average convergence within 
the Einstein radius is unity (see also section 3.1 of paper I). 
Suppose the model has stellar and dark matter components. 
Then the definition of the Einstein radius is 

Ktot(-Rein) = K.*(Rcin) + Kdm(-Rein) = 1 • (Bl) 

The stellar and dark matter mass fractions within the Ein- 
stein radius are, by definition, 

f'(T) \ — ^*(-^ cin ) fl (n \ _ ^dm(-Rcin) Con\ 

Ktot (,-flcin J Ktot (.-n-cinj 



Thus, equation (IBlfl can be rewritten as 

/:(i?cin)+/ dm ( J Rcin) = l (B3) 

On the one hand, this is entirely unsurprising: if the galaxy 
has stars and dark matter, the stellar and dark matter mass 
fractions must add up to unity. But on the other hand, 
this equation does illuminate our results in Fig. [5] Con- 
sider changing the dark matter component while keeping 
the stellar component fixed. As we vary the dark matter 
component, the Einstein radius will change. That in turn 
changes the first term in equation (|B3|) . which necessarily 
changes the second term, too. However, if we hold the stel- 
lar component fixed, the first term is a function of a single 
variable, namely the Einstein radius. Since / dm = 1 — fl 
that means the dark matter fraction can likewise depend 
only on a single variable, i? c i n . In other words, if we hold 
the stellar component fixed, then no matter what we do to 
the dark matter component, the changes to / dm propagate 
only through the Einstein radius, and we have a perfect one- 
to-one correspondence between i? c i n and /d m - 



